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ABSTRACT 


In the work reported herein, a simplified, uncoupled, zonal procedure is utilized to assess 
the capability of numerically simulating icing effects on a Boeing 727-200 aircraft. The 
computational approach combines potential flow plus boundary layer simulations by VSAERO 
for die un-iced aircraft forces and moments with Navier-Stokes simulations by NPARC for the 
incremental forces and moments due to iced components. These are compared with wind tunnel 
force and moment data, supplied by the Boeing Company, examining longitudinal flight 
characteristics. Grid refinement improved the local flow features over previously reported work 
with no appreciable difference in the incremental ice effect. The computed lift curve slope with 
and without empennage ice matches the experimental value to within 1 % , and the zero lift angle 
agrees to within 0.2 of a degree. The computed slope of the uniced and iced aircraft 
longitudinal stability curve is within about 2% of the test data. This work demonstrates the 
feasibility of a zonal method for the icing analysis of complete aircraft or isolated components 
within the linear angle of attack range. In fact, this zonal technique has allowed for the viscous 
analysis of a complete aircraft with ice which is currently not otherwise considered tractable. 
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1.0 INTRODUCTION 


The flow fields associated with complete aircraft configurations are extremely com- 
plicated and remain a challenge for applied computational aerodynamics. Linear methods, i.e., 
pure potential flow solutions that may include nonlinear wake shape effects, have become 
commonly used techniques for aerodynamic predictions and analysis; however, research 
scientists are now developing and applying nonlinear methods, thereby incorporating more 
accurate mathematical models of the actual flow physics, such as compressibility and viscosity. 
Although much progress has been made in solving full potential, Euler, and various forms of 
the approximate Navier-Stokes equations (NS), a zonal method of analysis is generally projected 
for the entire aircraft flow-field calculation because of the wide variation of physical scales of 
characteristic fluid phenomena. This method is necessary not only to achieve a closed 
calculation on the next generation of computers, but also to obtain the required numerical 
accuracy over large regions that must include fine meshes to resolve the viscous layers as well 
as large flow gradients. Computational requirements for simulating the effects of ice 
contamination are even more severe since the method must predict aircraft performance and 
stability degradation caused by the local formation of complex ice shapes. Furthermore, the 
method should also be practical and fast enough to be utilized in preliminary design studies to 
ascertain component sensitivity to icing. 

In this report, the work described involves a simplified, uncoupled zonal procedure which 
has been utilized to assess the capability of individual methods for simulating icing effects on 
a Boeing 727-200 aircraft. The zonal procedure includes the combined potential flow/boundary 
layer flow method (Program VSAERO 1 ) for the un-iced aircraft simulation, while a 
Navier-Stokes method (Program jnPARC 2 ) is utilized for the iced components of the aircraft. 
Wind tunnel data 3 used in this study for the un-iced and iced configurations were provided by 
tiie Boeing Company. The calculations described here include comparison with experiment for 
the longitudinal forces and moments. The grid resolution utilized in this study represents a 
refinement over previously reported work. Further results including lateral force and moment 
calculations may be found in an earlier report. 4 

The general method is described in Section 2.0, the geometric considerations are 
presented in Section 3.0, the results are detailed in Section 4.0, and the conclusions and 
recommendations are described in Section 5.0. A description of the routines added to NPARC 
along with a printout is included in Appendix A, listing of the utility program, PRPAMI, 
developed under this contract is included in Appendix B, a sample NPARC input file may be 
found in Appendix C, and VSAERO and NPARC manual updates are in Appendix D. The 
contracted work to enhance the panel model of the DHC-6-200 Twin Otter for fiiture tail-plane 
ice calculations is included in Appendix E. 
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2.0 METHODOLOGY 


Unlike the previous fully-coupled zonal procedure which provided for an iterative 
step, 5,6,7 the zonal aspect of this method is simplified to include an initial outer boundary flow 
condition only, which is provided by the potential flow code and utilized by the viscous flow 
component. At the viscous domain outer boundary (Fig. 1), velocities and pressure are 
specified and density is calculated from the pressure and free-stream sta gnatio n enthalpy. Since 
velocities, pressure and density are known, energy at the boundary can be determined. The 
viscous flow code uses these initial conditions in combination with the characteristic method to 
treat the outer surface boundary conditions. Further details pertaining to the implementation of 
flow boundary conditions may be found in the excellent discussion by Raj, et. al.*. 

The initial conditions generally involve effects from other aircraft features that are not 
included in the viscous domain. These may include, for example, lifting surface wake influence 
and fuselage blockage or upwash effects as well as the local surface boundary layer profile. In 
this manner, complications arising from an attempt to include other geometric features in the 
grid generation process are avoided. The final iced aircraft force and moment characteristics are 
estimated by superposition. The NPARC predicted incremental force and moment changes due 
to ice on individual aircraft components are combined with the VSAERO data for the un-iced 
aircraft. This process yields complete aircraft force and moment quantities under iced conditions. 
A zonal technique of this nature can allow for the viscous analysis of a complete aircraft with 
ice which is currently not otherwise considered tractable. 

2.1 VSAERO Technical Description 

VSAERO solves the Neumann problem of potential flow by converting T .aplace s partial 
differential equation into a boundary integral equation, as described by Morino 9 . The body and 
wake (i.e., boundary of the fluid) are then broken into finite size elements, or panels, each w ith 
an unknown potential value. Whereas, the body shape is known, the wake shape— determined 
by streamline trajectories from the downstream separation lines on the body— is not. An 
informed guess for the wake shape normally generates a useful first solution upon which 
VSAERO can improve iteratively. The displacement of the in viscid fluid boundary from solid 
surfaces caused by viscosity (namely, the boundary layer) can also be included by an iterative 
process. 

The panel method leads to a matrix equation of the order of the total number of surface 
panels. The matrix is dense and non-symmetric. The main features of the technique include 
generating the coefficients of the matrix, solving the matrix, and determining new wake shapes 
and viscous displacements. 
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2.2 NPARC Technical Description 

NPARC Version 1.0 (formerly PARC3D), is a multi-block code developed by Sverdrup 
Technology, Inc. 2 , and is based upon ARC3D 10 which was originally developed by Pulliam and 
Steger at NASA Ames Research Center. The method is based on the implicit approximate 
factorization algorithm of Beam and Warming 11 for unsteady compressible flow. For turbulent 
flow, the well-known Baldwin-Lomax model 12 is used for turbulence closure. This turbulence 
model also has been quite successful in the calculation of transonic flow, unless there is massive 
separation in the stream wise direction. 

NPARC contains several improvements over the original ARC3D code, the most notable 
being the multi-block capability. Also, the user has the option of specifying flow boundary 
conditions by means of namelist input. Generally, the modifications are directed toward 
improvement in user productivity. 

Several routines were modified or added to NPARC to accomplish the zonal component 
of this work. Generally the changes involved input/output of pertinent flow-field data and 
proper treatment of the far-field boundary condition. The outer domain zonal boundary 
condition implemented in this work was patterned after the freestream option available in 
NPARC (subroutine BCFAR, boundary condition type: 7) 2 . As previously mentioned, this 
boundary condition implementation is based upon the method of characteristics. The BCFAR1 
routine (derived from BCFAR), was developed to take advantage of the VSAERO flow 
information. Also, the contiguous block interface boundary condition (subroutine INTER70, 
boundary condition type: 70) was modified to remove the restriction that the computational 
coordinates must be increasing in the same physical direction for two overlapping grid blocks. 2 
The topological features of the t-tail grid-block structure did not conform to this requirement. 
A listing of the specific routines added to NPARC including a brief description are presented 
in Appendix A. 

2.3 Computational Procedure 

The step-by-step zonal procedure is outlined in Fig. 2. The user must first provide the 
basic model geometry in the form required by the particular grid generation package being used. 
Generally this process may be accomplished via a CAD database. The simple geometric tools 
available in VSAERO may be used for basic point re-distribution but in some cases the bi- 
quadratic interpolation has proven inadequate. Step two involves preparation of the volume grid 
within the flow domain of interest. In this work the grid generation process utilized the 
GRIDGEN system which results in the generation of a volume-grid data file in binary form. 

For projects of this computational magnitude, the hardware utilization generally follows 
that described in Fig. 3. The user conducts certain dements of the work on a local workstation 
and via network access, is able to perform the more computationally demanding tasks on a super 
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TASK 


SOFTWARE 


OUTPUT 


(1) PREPARE SURFACE 
GRID 

CAD AND/OR 
VSAERO 

VSAERO AND GRIDGEN 
SURFACE GEOMETRY 
DATA FILES 

(2) PREPARE VOLUME 
GRID 

GRIDGEN 

BINARY VOLUME 
GRID FILE 

(3) GENERATE NPARC 
RESTART DATA FILE 

PRPAMI 

NPARC RESTART DATA 
FILE 

(4) PREPARE VSAERO 
SCAN INPUT DATA 

PRPAMI 

VELOCITY SCAN INPUT 
DATA TO BE APPENDED 
TO VSAERO INPUT FILE 

(5) GENERATE OUTER 
FLOW B.C.'S AND 
UN-ICED CL, CM 

VSAERO 

VELOCITY SCAN OUTPUT 
DATA FILE 

(6) GENERATE VISCOUS 
FLOW SOLUTION 
INCLUDING CL, CM 

NPARC 

(AMI VERSION) 

NPARC RESTART DATA 
FILE 

(7) PREPARE GRAPHICS 
POST-PROCESSING 
DATA FILE 

PRPAMI 

0MNI3D GRAPHICS DATA 
FILE 


Fig. 2. Task S ummar y oi Zonal Ice Analysis Method. 
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HARDWARE UTILIZATION: 
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Fig. 3. Schematic of Hardware Utilization. 




computer. For investigations requiring the memory and computational speed of a main frame 
super computer, the transfer of the volume grid file should proceed in binary form. This will 
help preserve the numerical accuracy of the model inherent in the binary format. Subsequent use 
of the binary data on the super computer (usually 64 bit), may require translation from 32 bit 
to 64 bit format which may be accomplished with the following command: 


assign -F -f77 -N ieee ’filename’ 


This specific command is for a Cray computer running Umcos and should be used prior to each 
access of the binary data. 

The next task is normally conducted on the main frame computer and is involved with 
generation of the NPARC restart data file utilizing the utility program PRPAMI, which was 
developed under this contract. A copy of the code together with a brief description may be found 
in Appendix B. The GREDGEN volume grid file is required at this stage. The entire flow 
domain is initialized with frees tream flow conditions which are provided to PRPAMI 
interactively. 

The fourth step is the preparation of the VSAERO velocity scan input data. Again, 
PRPAMI is used in conjunction with the NPARC input file (INON5), and the restart file 
(REST. OLD). The user is responsible for specifying the flow boundary conditions in INON5, 
including the Type 8 zonal boundary condition developed for this work. (An example INON5 
data input file may be found in Appendix C). The PRPAMI code will assemble the grid point 
coordinates associated with all of the Type 8 regions and create a data file called SCAN.DAT 
in a format suitable for the VSAERO velocity input format. The user has the option of 
appending this data to the end of the VSAERO input files or using a separate data file for the 
velocity scan Mold 9 input. Modifications to the VSAERO and NPARC user’s manuals 
pertaining to this task may be found in Appendix D. 

The fifth task outlined in Fig. 2 involves the analysis of the complete model with 
VSAERO. The recommended approach is to conduct an incremental study of the model starting 
with rigid, specified wake trajectories and potential flow only. After verifying the solution 
integrity, the wakes may be relaxed and the viscous correction with the streamline-based 
boundary layer may then be added. This process will eventually provide the force and moment 
coefficient of the un-iced configuration of interest and yield the off-body velocity-scan flow 
calculation based upon the SCAN.DAT data file. The local surface boundary layer profile may 
be included as part of the NPARC flow ini tializa tion by setting the variable MEET to -3 in the 
VSAERO input file, (refer to Appendix D). 

The viscous flow solution for the isolated un-iced and iced model component is then 
determined as outlined under task six. As previously mentioned, the multi-block Reynold’s- 
averaged NS method NPARC was utilized for this purpose. Depending upon the particular 
application, the user may have to pay especially close attention to me convergence characteristics 
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of the NPARC solution procedure. This topic is treated in more detail in Section 3.0. A 
procedure to compute the force and moment coefficients has been added to NPARC as a 
requirement of this work. Based upon the user-specified boundary conditions, this module 
determines the force and moment quantities associated with each solid surface region in both 
body and wind axis systems. Total aerodynamic coefficients are included at the end of the region 
by region force and moment summary written to the NPARC output data file. A listing of the 
routines that comprise this force and moment module along with a brief description is included 
in Appendix A. 

Task seven should actually be utilized throughout the NPARC flow solution process 
initiated in task 6 to verify proper specification of flow boundary conditions and solution 
convergence characteristics. This step involves the generation of an OMNDD compatible 
graphics data file with the utility code PRPAMI. The end result is a formatted data file which 
can be transferred to the OMNDD platform for graphical processing. 
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3.0 GEOMETRY MODELING 


This section details the work involved in the preparation of the Boeing 727-200 data case 
which was selected for validation of the proposed zonal concept. The current effort involves the 
analysis of the model at multiple angles of attack and zero sideslip in an effort to estimate the 
aerodynamic degradation due to the presence of leading-edge ice on the horizontal and vertical 
tail. 


3. 1 Panel Model 

A VSAERO model of a complete Boeing 727-200 was prepared for this study. The 
wing was constructed with a jig twist distribution, which was consistent with the Boeing wind 
tunnel model (Fig. 1). The horizontal tail setting was -4 degrees with respect to the fuselage 
reference line. The center engine inlet was modified by removing the ducting inside the 
highlight and adding an elliptical faring from the highlight forward. The outboard engine was 
changed to a flow-through nacelle by extending the inlet duct aft to the trailing edge of the 
nacelle. Unfortunately, Boeing was not able to provide pictures of the support sting that was 
attached to the model at the aft end of the fuselage. In lieu of this information, a wake was 
attached to the trailing edge of the center engine exhaust to simulate the effects of the sting; 
therefore the accuracy of the modeling of the sting compared with the actual experiment is 
unknown at this time. 

The VSAERO panelization of the starboard side of the aircraft included 4,105 body and 
3,899 wake panels in order to simulate the symmetric flight conditions described in this report. 

3.2 Grid Generation 

Structured multi-block grids were developed individually about the isolated t-tail for both 
iced and un-iced flight conditions. Because of the zonal approach employed here, the NS domain 
for both grids was bounded by the region described in Fig. 1. The outer domain boundary 
shown here is about one semi-span from the horizontal tail. Further, the lower domain boundary 
includes a section of the vertical tail above the center engine. As previously described, the flow 
information on the outer faces of these volumes are provided by the VSAERO code including 
an initial estimation of the boundary layer growth on the vertical tail section. 

The minimum grid spacing normal to the body surface is consistent with a y+ of order 
10 (minimum AZ/C = 0.0001). Due to the presence of the vertical tail, an H-0 grid topology 
was used to discretize the flow field (Fig. 4). This representation generally does not resolve the 
leading-edge region as well as a C-type topology but the t-tail geometric constraints limit the 
topological choices. The three-block grids about the un-iced and iced t-tail consisted of over 
690,000 grid points (194x73x49) and 815,000 (228x73x49) grid points, respectively (See 
Appendix C). This represents an increase of nearly 90,000 grid points over the grid resolution 
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Schematic of T-tail H-0 Grid Topology 
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Fig. 4. Schematic of T-tail H-0 Grid Topology. 



used in previously reported work 413 . The region in the vicinity of the horizontal tail leading-edge 
has been refined in the work presented here as shown in Fig. 5. Geometry details of the 
leading-edge ice shape for the horizontal tail are presented in Section 4.0. 

Volume grids for both the un-iced and iced t-tail configurations were generated using 
GREDGEN3D Version 8.5 14 . Special modifications were made by MDA Engineering, Inc. to 
improve the control of orthogonality on selected faces. Specifically, the control functions were 
averaged across the leading edge of the horizontal tail (£ direction), the leading edge of the 
vertical tail (£ direction), and the tip of the vertical tail (f direction). Also, the foreground 
control functions were only calculated on the region of the tip block that corresponded to the tip 
surface of the horizontal tail. It is anticipated that these features will be included in version 9 
of GREDGEN3D 15 . 

The geometry manipulation required in this work was continually plagued with questions 
regarding accurate representation of the basic surface. The initial surface definition, based upon 
tabular data supplied by Boeing, as well as the single precision software components of the 
surface and grid generation processes, contributed to these difficulties. Clearly, grids that 
achieve y + values of 0(1) are desirable for the viscous simulations reported here; however, this 
level of exactness requires surface definition to a minimum of eight places of accuracy and 
software tools capable of double precision computations. 
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ri 6 . 5. Comparison of Original and Revised Grid in the 
Vicinity of the Horizontal Tail Tip Section. 
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4.0 RESULTS 


A 0.075-scale model of the Boeing 727-200 was tested in January of 1981 in the Boeing 
VERTOL Tunnel. Wind tunnel data was obtained at a Mach number of approximately 0. 12 and 
a Reynolds number of 1.24 million, based on wing reference chord. Longitudinal forces and 
moments were measured for the un-iced cruise configuration and for various combinations of 
iced leading edges on the horizontal and vertical tail. A summary of the specific flight conditions 
reported here is included in Table 1 below. This section presents aircraft force and moment 
comparisons to experiment which have been previously reported 413 , as well as detailed flow 
features for the iced t-tail at zero degrees angle of attack. 


CASE 

DESCRIPTION 

CODE 

ALPHA 

(deg) 

YAW 

(deg) 

MACH 

Re 

Un-iced Aircraft 

VSAERO 

-4, 0,2,4 

0.0 

0.12 

1.2 x E06 

Un-iced T-Tail 

NPARC 

0,4 

0.0 

0.12 

1.2 x E06 

Iced T-Tail 

NPARC 

0,4 

0.0 

0.12 

1.2 x E06 


Table 1. Summary of Flow Calculations. 


In previous work 16 , it was not possible to generate an acceptable t-tail grid which would 
meet the single-block assumption and thin layer approximation of ARC3D. Conse-quently, 
longitudinal icing calculations included ice on the horizontal tail only. Icing effects, created by 
ice on the vertical tail and ice on the outboard Kruegers, were included in the test data, but were 
neglected in the longitudinal force and moment calculations. 

In the current work, the multi-block capability of NPARC removes these earlier 
restrictions and allows for an improved representation of the empennage. The longitudinal 
characteristics are examined with leading-edge ice on both horizontal and vertical tail control 
surfaces as required. Based on practical considerations, the ice on the Krueger flap has been 
ignored in the present calculations. 
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4. 1 Un-Iced Aircraft Simulation 


Since the Mach number in the test was essentially negligible, all VSAERO runs were 
made without compressibility corrections. Longitudinal calculations included angle-of-attack 
settings of -4, 0, 2 and 4 degrees. Computed and measured force coefficients were reduced with 
a reference chord of 180.7 inches, a reference area of 244,800 square inches, and a half span 
of 648 inches. Moments were referenced to fuselage station 905.3 and waterline 184.2 inches. 

The viscous correction for the un-iced VSAERO analysis was accounted for through the 
coupled, streamline-based integral boundary layer procedure. In the boundary layer 
calculations, transition was set to occur as soon as the local Reynolds number, based upon 
momentum thickness, reached a value of 200. Without other test information regarding 
transition behavior, this method has been found to provide an appropriate boundary layer 
simulation for such wind tunnel test conditions. Three viscid-in viscid correction iterations were 
computed for each flight condition. Computed contours of constant pressure coefficient are 
presented in Fig. 1 for zero degrees angle of attack. 

Figs. 6 and 7 compare the computed lift and pitching moment characteristics with the 
experimental data. The VSAERO inviscid computation is included hoe for q ualit ative 
verification of the boundary layer effects. Within the linear range, the VSAERO boundary layer 
si mu la t i o n matches the slope of the experimentally measured lift curve to within 1 % and the zero 
lift angle agrees to within about one-tenth of a degree. This comparison is greatly improved 
over that previously reported in Ref. 16, which exhibited an error of approximately three-fourths 
of a degree. After further consultation with Boeing, it was determined that the defined wing 
twist distribution utilized in the previous work did not match the model specification; 
consequently, the current results represent the correct wing twist distribution. The calculated 
pitching moment slope of -0.032 compares well to the test value of -0.029 per degree and the 
computed zero pitching moment angle of attack is within 0.33 degree of the test value. 
Discussion of the un-iced aircraft longitudinal stability curve is included in the next section. 

4.2 Iced Aircraft Simulation 

In order to compute the force and moment increments of the aircraft due to ice, the 
viscous flow characteristics of the Boeing 727-200 empennage with and without leading-edge 
glaze ice was examined. Table 1 provides a detailed summary of the NPARC computations. 
Each of the NPARC computations required approximately 2.4xE-0.05 sec. /grid point/iteration 
and 23 Mw of memory on a CRAY-YMP Model-E. 
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Boeing 727-200 Lift Comporison 
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Fig. 6. Comparison of Measured to Calculated Aircraft Lift. 






Boeing 727-200 Pitching Moment Comparison 
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Fig. 7. Comparison of Measured to Calculated Aircraft Pitching Moment. 



In all of the NPARC simulations, the flow was assumed to be fully turbulent at a model 
aircraft Reynolds number of 1.24 million. In these simulations, NPARC required from 10,000 
to 15,000 iterations to achieve residual reductions on the order of 2 to 3 for a CFL number of 
2 (see for example, Fig. 8). In practice, residual reduction of 3 to 4 orders is the generally 
accepted convergence criterion. In this study however, various factors have contributed to the 
poor residual convergence behavior. First, the fine grid resolution tends to restrict the maximum 
CFL number or time-step size which results in an increased number of iterations. 17 Second, the 
multi-block nature of the domain decomposition adversely affects solution convergence due to 
the time lag of the boundary information. 1 * Third, the low subsonic Mach number is at or near 
the operational limit of the compressible flow equation set. 17 Further testing is required to 
quantify the relative contribution of each factor and to suggest means by which the residual 
convergence can be accelerated. 

In addition to residual history, the integrated lift and pitching moment coefficients were 
monitored as an aid in determining the convergence level. The history of lift and pitching 
moment coefficient is shown for the iced and un-iced t-tail configuration in Fig. 9 for an angle 
of attack of zero degrees. The unsteady behavior previously noted 413 , is not obvious in the 
current calculations. It is apparent that the grid refinement and/or code differences between 
PARC3D and NPARC have contributed to this improvement in computed flow characteristics. 
As shown by the time histories of lift and pitching moment, the iced t-tail required more 
iterations to insure adequate convergence levels in comparison to the un-iced geometry. 

A comparison of the un-iced and iced horizontal tail mid-span sectional pressure 
distribution at 0° angle of attack is presented in Fig. 10. As indicated, the presence of 
leading-edge ice generates additional down force resulting in an increased aircraft nose-up 
pitching moment. A close examination of the iced airfoil indicates that the ice shape produces 
a negative droop of the leading-edge. Thus the ice shape modifies, locally, the airfoil camber 
to produce the additional loading. There is also an apparent contribution from the effect of 
leading-edge ice on the boundary layer which would require further study to identify. 

The horizontal tail reversed flow region in the vicinity of the upper and lower ice horn 
is presented in Fig. 11 for both original and current grid resolution. It is apparent that the flow 
solution on the refined grid is improved at least in a qualitative sense over the result on the 
original grid. A close-up view of the lower hom clarifies this point further (Fig. 12). The solid 
horizontal line in both figures marks the boundary between the upper and lower grid blocks. 
Note especially the flow vectors in the vicinity of the ice hom leading-edge and block interface. 
The refined grid result shows marked improvement over the original grid calculation. 

The computed increment in lift and moment coefficients due to the presence of ice at zero 
degrees angle of attack is -0.0025 and 0.0122 respectively; while the corresponding increments 
based upon the original grid density are -0.0017 and 0.0095. This does not appreciably modify 
the previously reported force and moment results. Furthermore, this result tends to alleviate the 
earlier concerns regarding repeatability due to the relatively small force and moment increments 
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Fig. 10. NPARC Computed Surface Pressure Distribution 
at Mid-Span of the Horizontal Tail. 
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Fig. 11(a). Flow Velocity Vectors in the Vicinity of the Iced Horizontal Tail 
Leading Edge for the Original Grid. Flow is from Right to Left. 
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:ity Vectors in the Vicinity of the Iced Horizontal Tail 
Ige for the Revised Grid. Flow is from Right to Left. 
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Fig. 12(a). Close-up View of Flow Velocity Vectors in the Vicinity of the Lower Lobe of 
Iced Horizontal Tail for the Original Grid. Flow is from Right to Left. 




due to ice 4,13 . In feet, the detailed force and moment comparisons presented in this report are 
based on calculations utilizing the original grid resolution. 


The computed aircraft lift and pitching moment characteristics are compared with experi- 
mental data in Figs. 13 and 14 and a quantitative breakdown is presented below in Tables 2 and 
3. For comparison purposes, the test data characteristics were estimated by restricting the region 
of interest to the linear angle-of-attack range. The computed lift curve slope with ice matches 
the experimental value to within 1 %, and the zero lift angle agrees to within two-tenths of a 
degree. In addition, the computed pitching moment curve slope of -0.035 correlates with the 
test value of -0.032 per degree. 


test 

CALCULATION 

UN-ICED 

ICED 

dC L /da 
(per deg) 

a. 

(deg) 

Cu 

X 


X 


.077 

-1.48 

.114 

X 



X 

. 077 

-1.44 

.111 


X 

X 


.084 

-1.36 

.114 


X 


X 

.085 

-1.29 

.110 


Table 2. Detailed Comparison of the Aircraft longitudinal Characteristics:C L vs, a. 
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Boeing 727-200 Lift Comparison 
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Fig. 13. Comparison of Measured to Calculated Impact of Ice on Aircraft Lift Coefficient. 
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Fig. 14. Comparison of Measured to Calculated Impact 
of Ice on Pitching Moment Coefficient. 



TEST 

CALCULATION 

UN-ICED 

ICED 

dC M y/da 
(per deg) 

ABSCISSA 

INTERCEPT 

(deg) 

X 


X 


-.029 

8.15 

X 



X 

-.032 

7.94 


X 

X 


-.032 

7.84 


X 


X 

-.035 

7.80 


Table 3. Detailed Comparison of the Aircraft Longitudinal Characteristics rC^ vs. a. 


As indicated, the computed absolute slope and intercepts compare well with the test data as well 
as the increments due to ice. In particular, several incremental quantities may be estimated from the 
tabulated data. The experimental increment due to ice in lift curve slope, (AC U ), zero lift angle of attack, 
(aj, and lift at zero angle of attack, (AC,J, of 0.000, 0.04 and -0.003 compares very well with their 
respective computed values of 0.001, 0.07 and -0.004 (see Table 2). The experimental increment due 
to ice in pitching moment curve slope, (AC^J, and abscissa intercept, (Aa), of -0.003 and -0.21 
correlates well with their respective computed values of -0.003 and -0.04 (see Table 3). Evidently, for 
the longitudinal case, the additional effects of the ice on the outboard Kruegers that was included in the 
experiment are indeed small for the examined angle of attack range. 

The aircraft longitudinal stability curve is graphically presented in Fig. 15 and a detailed 
comparison of the computed and test stability data are presented below in Table 4. As shown, the 
computed slopes for the un-iced and iced aircraft are within about 2% of the test data. Further, based 
on estimates derived from Table 4, the experimental increment due to ice in stability curve 
slope, (dCwy/dCO and pitching moment at zero lift (AC,*,), of 0.16 and -0.01 compares favorably with 
the computed values of 0. 18 and -0.003 respectively. In this case, the addition of ice to the vertical and 
horizontal tail leading-edge is a destabilizing influence. 
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Fig. 15. Comparison of Measured to Calculated Impact of Ice on Aircraft Stability. 





TEST 

CALCULATION 

UN-ICED 

ICED 

dCm*/d* 
(per deg) 

Cmt 

X 


X 


-2.63 

.761 

X 



X 

-2.47 

.751 


X 

X 


-2.60 

.773 


X 


X 

-2.42 

.770 


Table 4. Detailed Comparison of the Aircraft Longitudinal Characteristics: C^ vs. C L . 


Generally, computed increments due to ice are in the proper direction and are of the same 
magnitude as the test data for the aircraft longitudinal characteristics. The small relative magnitude of 
these increments, though, makes it difficult to discuss error in terms of percent. 
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5.0 CONCLUSIONS AND RECOMMENDATIONS 


In the work reported here, a simplified, uncoupled, zonal procedure is utilized to assess the 
capability of numerically simulating icing effects on a Boeing 727-200 aircraft. The computational 
approach combines potential flow plus boundary layer simulations by VSAERO for the un-iced aircraft 
forces and moments, with NS simulations by NPARC for the incremental forces and moments due to iced 
components. The computed lift curve slope with and without empennage ice matches the experimental 
value to within 1%, and the zero lift angle agrees to within two-tenths of a degree. Furthermore, the 
computed slope of the un-iced and iced aircraft longitudinal stability curve is within about 2% of the test 
data. These results demonstrate the feasibility of a zonal method for the icing analysis of complete aircraft 
or isolated components within the linear angle of attack range. In fact, this zonal technique has allowed 
for the viscous analysis of a complete aircraft with ice which is currently not otherwise considered 
tractable. 

Future work should include investigation of solution convergence acceleration methods in order 
to reduce the current CPU requirements, and double precision software in order to improve the geometric 
accuracy. Additionally, an experimental data base that includes surface definition to a minimum 0 f eight 
places of accuracy, and an incremental buildup of the iced components as well as surface pressure 
information should be available in order to thoroughly evaluate the CFD prediction capability. 
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APPENDIX A 


The following subroutines have been added to the NS code NPARC version 1.0 and were 
developed under NAS3-26310. These modifications pertain to the zonal aerodynamic analysis of aircraft 
or aircraft components with ice. 


Subroutine 

Name Description 


RFIXED 

BCFAR1 

SETVEL 

FORMOM 

PATCH 

CPCALC 

CORNR 

ANUZ 


Opens outer boundary condition data file PARC.SCN as fortran unit 32. Loads 
uo, vo, wo, cpo arrays. Called from INlllA if type 8 boundary condition 
specified in input. 

Patterned after routine BCFAR. Sets zonal boundary condition based upon uo, 
vo, wo, cpo information. Called from routine BC. 

Sets local flow parameters from uo, vo, wo, cpo arrays. Called from BCFAR1 . 

Computes the force and moment coefficients for all type 60 and 61 regions. 
Called from main routine. 

Forms surface panel element information such as, area and normal vector. 
Called from FORMOM. 

Computes local surface panel element pressure coefficient. Called from 
FORMOM. 

Assembles flow variables at each comer of a surface panel element. Called from 
CPCALC. 

Performs loads summation over each surface region and prints out result. Called 
from FORMOM. 
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SUBROUTINE RFIXED 

C** ************************************************************* ******** 
cdan New routine added in December 1992 for VSAERO Vscan data, 
c*************************************************************** ******** 
include ' npami . inc ' 

PARAMETER (NSCNMX=JMX*KMX*IMX) 

COMMON/ SCAN/UO (NBX , JMX, KMX , IMX) , VO(NBX, JMX,KMX,EMX) , 

& WO (NBX, JMX, KMX, LMX) ,CPO (NBX, JMX, KMX, IHX) 

OPEN (UNIT™ 3 2 , FILE™' PARC. SCN' , FORM™ ' FORMATTED' ,ERR=802, 

& STATUS™ ' OLD ' ) 

REWIND 32 
NPT=»0 

C SET DEFAULTS 

DO 4 N=1 , NBX 
DO 4 J=1 , JMX 
DO 4 K=1 , KMX 
DO 4 L=1 , LMX 
UO(N, J , K, L) =*501. 

4 CONTINUE 

C READ FIRST TWO DUMMY RECORDS . . . 

READ(32 , *) 

READ(32 , *) 

DO 5 J=1 , NSCNMX 

READ(32,*,END=801,ERR=5)U,V,W,CP,NB,JP,KP,LP 

IF ( JP . EQ . 0 ) GO TO 5 

NPT=NPT+1 

UO(NB,JP,KP,LP)™U‘ 

VO(NB,JP,KP,LP)«V 
WO (NB , JP, KP , LP) «W 
CPO (NB , JP , KP , LP) -CP 

5 CONTINUE 

WRITE (6,*) ' ERROR IN RFIXED, ALL BC POINTS HAVE NOT BEEN READ' 
WRITE (6,*)' INCREASE NSCNMX PARAMETER' 

STOP 

801 CONTINUE 

WRITE (*, *) '*** TOTAL NUMBER OF FIXED BC POINTS READ. . .',NPT 
RETURN 

802 WRITE ( * , * ) ' *** ERROR OPENING UNIT 32 ***' 

STOP 

END 

SUBROUTINE BCFAR1 (ID, J,KA,KB,LA,LB, ISGN, PINF , TINF , 

* XX,XY,X2, YX,YY,YZ, ZX,ZY,ZZ, Q ) 

C**** ********************************************** ********************* 
cdan New routine added in December 1992 for VSAERO Vscan data. 

C*********************************************************************** 

include 'npami. inc' 

COMMON/INTEG/ IBF(NM+1) , FNMACH(NM) , f ill (58*NM-1+12*NDIF) 
COMMON/VIS FAC/ ALPHA, BETA , XMACH , C2B, C2BP, VRAT,RE 
COMMON/WAYOUT/ AO , RO , UO , VO , WO , EO , PO 
COMMON/MB LOCK/NB LOCK , MB , MBNEXT 
DIMENSION Q ( JMAX, KMAX,IMAX, 6) , 

* XX(Jmax,Rmax,Lmax) , 

* XY(Jmax,Kmax,Lmax) , 

* XZ(Jmax,Kmax,Lmax) , 

* YX(Jmax,Kmax,Lmax) , 

* YY(Jmax,Kmax,Lmax) , 


36 



* YZ ( Jmax, Kmax, Lznax) , 

* ZX (Jmax , Kmax, Lmax) , 

* ZY (Jmax, Kmax, Lmax) , 

* ZZ (Jmax, Kmax, Lmax) 

s 

C OUTER BOUNDARY FLOW CONDITIONS 

C CAUTION: ONLY SATISFACTORY FOR INCOMPRESSIBLE FLOW CONDITIONS 
AO = SQRT( TINF ) 

RO = GAMMA* PINF/TINF 


JP=J+ISGN 
SGN=ISGN 
DO 310 L=LA,LB 
DO 301 1=1 , NM+1 
IBF ( I ) =0 
301 CONTINUE 
C 

IF (ID.EQ.l) THEN 
DO 303 K=KA, KB 

ZNORM = SGN * SQRT ( ( (XX (J,K,L)**2) +XY ( J,K, L) **2+XZ (J,K,L)**2) ) 
DCX = XX(J,K,L) /ZNORM 

DCY = XY(J,K,L) /ZNORM 

DCZ = XZ(J,K,L) /ZNORM 

C EXTRAPOLATE TO BOUNDARY 

RI = Q ( JP, K, L, 1) *Q ( JP,K, L, 6) 

UI = Q(JP,K,L,2)/Q(JP,K,L,1) 

VI = Q(JP,K, L, 3)/Q ( JP,K, L, 1) 

WI = Q(JP,K,L,4)/Q(JP,K,L,1) 

El = Q ( JP, K, L, 5) *Q ( JP,K, L, 6) 

PI = GAMI* ( El -0. 5*RI* (UI**2 + VI**2 + WI**2) ) 

UNI = DCX*UI +DCY*VI +DCZ*WI 

AI = SQRT ( GAMMA*PI/RI ) 

C OUTER BOUNDARY NORMAL VELOCITY 

CALL SETVEL(MB, ID, J,K,L,U0 , VO ,W0 ,CP0) 

U0=U0*XMACH 
V0=V0*XMACH 
W0=W0*XMACH 
P0=PINF+ . 5 * CPO * XMACH 

EO = PINF*GM1R + . 5*R0* (U0**2 + V0**2 + W0**2) 

ctmp 

c write (6,*) '*** bcfarl: ' ,ID, j , jp,k, I,u0,v0,w0,p0 

UNO = DCX*U0 +DCY*V0 +DCZ*W0 
C COMPUTE SPEED OF SOUND AND NORMAL VELOCITY COMPONENT 
C FROM RIEMANN INVARIANTS 

Rl = UNO + ( 2 . 0*GM1R) *A0 

R2 = UNI - (2.0 *GM1R) *AI 

A = 0.25* (R1-R2) *GAMI 

UN = 0 . 50* (R1+R2) 

XMN = UN/A 

FNMACH(K) =XMN 

303 CONTINUE 

ELSE IF (ID.EQ.2) THEN 
DO 304 K=KA KB 

ZNORM = SGN * SQRT ( ( (YX (K, J, L) **2) +YY (K, J,L) **2+YZ (K, J,L) **2) ) 
DCX = YX(K,J,L) /ZNORM 

DCY = YY(K,J,L) /ZNORM 

DCZ = YZ(K,J,L) /ZNORM 
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C EXTRAPOLATE TO BOUNDARY 

RI = Q(K, JP,L, 1) *Q(K, JP,L, 6) 

UI = Q(K,JP,L,2)/Q(K, JP,L, 1) 

VI = Q (K, JP, L, 3) /Q(K, JP, L, 1) 

WI = Q (K, JP, L, 4) /Q(K, JP, L, 1) 

El « Q (K, JP, L, 5) *Q(K f JP, L, 6) 

PI = GAMI* ( El -0.5*RI*(UI**2 + VI**2 + WI**2) ) 

UNI = DCX*UI +DCY*VI +DCZ*WI 

AI = SQRT( GAMMA*PI/RI ) 

C OUTER BOUNDARY NORMAL VELOCITY 

CALL SETVEL(MB,ID, J, K, L,UO , VO , WO , CPO) 

UO=UO*XMACH 
VO=VO*XMACH 
WO=WO*XMACH 
PO=PINF+ . 5*CP0*XMACH 

EO = PINF*GM1R + . 5*R0* (U0**2 + V0**2 + W0**2) 

ctmp 

c write(6, *) '*** bcfarl: ' , ID, j , jp,k, l,uO,vO,wO,pO 

UNO = DCX*UO +DCY*VO +DCZ*WO 
C COMPUTE SPEED OF SOUND AND NORMAL VELOCITY COMPONENT 
C FROM RIEMANN INVARIANTS 

Rl = UNO + ( 2 . 0*GM1R) *AO 

R2 « UNI - (2 . 0*GM1R) *AI 

A = 0. 25* (Rl— R2 ) *GAMI 

UN = 0.50* (R1+R2 ) 

XMN = UN/A 

FNMACH(K) =XMN 

7 04 CONTINUE 

ELSE IF (ID.EQ.3) THEN 
DO 305 K=KA,KB 

ZNORM = SGN * SQRT( ( (ZX(K,L, J) **2)+ZY(K,L, J) **2+ZZ(K,L, J) **2) ) 
DCX = ZX(K,L,J) /ZNORM 

DCY = ZY(K,L,J) /ZNORM 

DCZ = ZZ(K,L,J) /ZNORM 

C EXTRAPOLATE TO BOUNDARY 

RI = Q(K,L, JP, 1) *Q(K,L, JP, 6) 

UI = Q(K,L,JP,2)/Q(K,L,JP,1) 

VI = Q(K,L,JP,3)/Q(K,L,JP,1) 

WI = Q(K,L,JP,4)/Q(K,L,JP,1) 

El = Q(K,L, JP, 5) *Q(K,L, JP, 6) 

PI => GAMI* ( El -0. 5*RI* (UI**2 + VI**2 + WI**2) ) 

UNI = DCX*UI + DCY* VI +DCZ*WI 

AI - SQRT ( GAMMA* PI/RI ) 

C OUTER BOUNDARY NORMAL VELOCITY 

CALL SETVEL(MB, ID, J,K,L,U0,V0,W0,CP0) 

UO“UO*XMACH 
VO=VO*XMACH 
W0-W0*XMACH 
PO«=PINF+ . 5*CP0*XMACH 

EO = PINF*GM1R + . 5*R0* (U0**2 + V0**2 + W0**2) 

ctmp 

c write(6,*) '*** bcfarl: ', ID, j , jp,k,l,uO,vO,wO,pO 

UNO = DCX*U0 +DCY*V0 +DCZ*W0 

COMPUTE SPEED OF SOUND AND NORMAL VELOCITY COMPONENT 
C FROM RIEMANN INVARIANTS 

Rl = UNO + ( 2 . 0*GM1R) *A0 
R2 = UNI - (2 . 0*GM1R) *AI 

A = 0. 25* (R1-R2) *GAMI 
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UN = 0 . 50* (R1+R2 ) 

XMN = UN/A 


FNMACH (K) =XMN 

305 CONTINUE 
END IF 

DO 307 K=KA, KB 
IF (FNMACH (K) . GT.1.0) THEN 
IBF (K) =1 ‘SUPERSONIC INFLOW 

ELSE IF (FNMACH (K) .GT.0.0) THEN 
IBF (K) =2 ! SUBSONIC INFLOW 

ELSE IF (FNMACH(K) .GT.-1.0) THEN 
IBF (K) =3 ! SUBSONIC OUTFLOW 

ELSE 

IBF (K) =4 [SUPERSONIC OUTFLOW 

END IF 
307 CONTINUE 
Kl—KA 

IBFF=IBF (KA) 

DO 309 K=KA+1,KB+1 
IF (IBF(K) .NE.IBFF) THEN 
K2=K-1 

IF (IBFF.EQ. 1) THEN 

CALL FARFIX (ID, J, K1 ,K2 , L, BIGA(IPWR) ) 

ELSE IF (IBFF.EQ. 2) THEN 
CALL FARIN ( ID, J , JP, K1 , K2 , L,SGN, 

* BIGA(IPWXX) , BIGA(IPWXY) , BIGA(IPWXZ) , 

* BIGA(IPWYX) , BIGA (IPWYY) ,BIGA(IFWYZ) , 

* BIGA(IPWZX) , BIGA ( IPWZ Y ) ,BIGA(IFWZZ) , 

* BIGA (IPWR ) ) 

ELSE IF (IBFF.EQ. 3) THEN 

CALL FAROUT (ID, J , JP, K1 , K2 , L,SGN, 

* BIGA (IPWXX) , BIGA(IPWXY) ,BIGA(IPWXZ) , 

* BIGA(IPWYX) , BIGA (IPWYY) ,BIGA(IPWYZ) r 

* BIGA(IPWZX) , BIGA ( IPWZ Y) , BIGA ( IPWZ Z) , 

* BIGA (IPWR ) ) 

ELSE IF (IBFF.EQ. 4) THEN 

CALL OUTSUP ( ID , J , JP , K1 , K2 , L , BIGA(IPWR) ) 

END IF 
K1~K 

IBFF=IBF(K) 

END IF 

309 CONTINUE 

310 CONTINUE 
END 


SUBROUTINE SETVEL (MB , ID , JPP , K , L , UO , VO , WO , CPO ) 

C** ***************** **************************************************** 

cdan New routine added in December 1992 for VSAERO Vscan data. 


**** 


include ' npami . inc 9 

COMMON/ SCAN/UO (NBX, JMX,KMX, UfX) , VO (NBX, JMX, KMX, 123" , 
& WO (NBX, JMX, KMX, LMX) , CPO (NBX , JMX, KMX, LMX) 


U0=1001. 

IF (ID. EQ. 1) THEN 
UO=UO (MB, JPP, K, L) 
VO=VO (MB, JPP, K, L) 
WO=WO(MB,JPP,K,L) 
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o o 


CPO-CPO (MB , JPP , K , L) 

ELSEIF ( ID . EQ . 2 ) THEN 
UO=UO (MB , K, JPP, L) 

VO=VO (MB, K, JPP, L) 

WO*WO (MB , K , JPP , L) 

CPO=CPO (MB , K , JPP , L) 

ELSEIF (ID. EQ . 3 ) THEN 
UO=UO(MB,K, L, JPP) 

VO=VO (MB , K , L , JPP) 

WO=WO ( MB , K , L , JPP) 

CPO=CPO (MB , K , L, JPP) 

ENDIF 

IF(UO.GT. 1000. ) THEN 

WRITE(6, *) '*** ERROR. . . VELOCITY NOT SET ***' 

WRITE (6,*)' INDICES DO NOT MATCH SCAN DATA' 

WRITE ( 6 , * ) ' MB, ID, JPP, K, L, ' ,MB, ID, JPP,K,L 
STOP 

ELSEIF (DO . GT. 500 . ) THEN 

WRITE(6,*) '*** ERROR... VELOCITY NOT SET ***' 

WRITE (6,*) ' SCAN VELOCITY NOT AVAILABLE' 

WRITE (6 ,*) * MB, ID, JPP, K, L, ' ,MB, ID, JPP,K,L 
STOP 
ELSE 
RETURN 
ENDIF 
END 

SUBROUTINE FORMOM 

C ***** ***************** ********* *************** ****************** ******* 

'"\n Routine added for force & moment calc. 7/93 

.■A******************************************************************** 

C 

include 'npami.inc' 

PARAMETER (MAXB=*4 *NIP) 

COMMON/MBLOCK/NBLOCK , MB , MB NEXT 
COMMON/TIME/ DT, IVARDT 

COMMON/FORCE/ CBAR , SSPAN , SREF , REFMX , REFMY , REFMZ 
COMMON/VISFAC/ ALPHA, BETA, XMACH, C2B, C2BP, VRAT,RE 
COMMON/KEPS / NRLX, NTURB, ORDER, SIGK, SIGEPS , CK1, CK2 , CK3 , CK4 , CK5 
COMMON/INTEG/JL (MAXB) , KL(MAXB) , LL(MAXB) 

DIMENSION Q(JMX,KMX,LMX,5) , X ( JMX, KMX, LMX) , 

1 Y( JMX, KMX, LMX) , Z (JMX, KMX, MX) ,XC(NIP) ,YC(NIP) ,ZC(NIP) , 

2 CPP(NIP) , AREA (NIP) ,ANT(NIP,3) 

NAMELIST/BOUNDS/ NJSEG , JUNE , JKLOW , JKHIGH , JLLOW , JLHIGH , 

' NKSEG, KLINE, KJLOW, KJHIGH,KLLOW,KLHIGH, 

' NLSEG , LLINE , LJLOW , L7HIGH , LKLOW , LKHIGH , 

» JTYPE , JSIGN , PRESS J , TEMPJ , INTERJ , 

> KTYPE , KSIGN , PRESSK, TEMPK, INTERK , 

r LTYPE , LSIGN , PRESSL, TEMPL, INTERL , 

' JEDGE , KEDGE , LEDGE , JDIR , KDIR, LDIR, 

» NORMJ , NORMK , NORML 

NTI = 5 
REWIND NTI 
NBPT - 0 
NPAT = 0 

WRITE (6, 800) 

WRITE (6, 801) CBAR, SREF, SSPAN, REFMX, REFMY, REFMZ, XMACH, ALPHA, BETA 
RETRIEVE THE 3D SOLUTION DATA 
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REWIND 4 

READ (4 )NC, GAMMA 

DO 50 NB=1 , NBLOCK 

READ(4) JL(NB) ,KL(NB) , LL(NB) 

READ (4) 

* ( ( (X ( J, K, L) , 

* J=l, JL(NB) ) , K=1 , KL(NB) ) ,L=1,LL(NB) ) , 

* (((Y(J,K,L), 

* J=1,JL(NB) ) , K=1 , KL(NB) ) ,L=1,LL(NB) ) , 

* ( ( (Z (J, K, L) , 

* J=l, JL(NB) ) , K=1 , KL(NB) ) , L=1 , LL(NB) ) 

READ(4 ) ( ( ( (Q(J,K,L,N) , J=1,JL(NB)) ,K=1,KL(NB)) , 

* L=1 , LL(NB) ) , N=1 , 5) 

IF(NC.ge.NTORB) 

* READ(4) (BIGA(J) ,J=IPWAK , IPWAK +NXYZ-1) , 

* (BIGA(J) , J=IPWEPS , IPWEPS+NXYZ-1) , 

* (BIGA(J) , J=IPWTMU, IPWTMU+NXYZ— 1) 

READ (NTI, BOUNDS) I Cray version of NAMELIST 

c READ ( NTI, NML= BOUNDS) ! IRIS version of NAMELIST 

DO 26 J=l, NJSEG 

IF(JTYPE(J) .EQ.60.0R. JTYPE(J) .EQ. 61) THEN 
IORDR=l 
NPAT=NPAT+ 1 

CALL PATCH (NORMJ(J) , IORDR, JKLOW(J) , JKHIGH(J) ,JLLOW(J) , 

1 JLHIGH(J) , JLINE ( J) , X,Y, Z , AREA, ANT, XC,YC, ZC) 

CALL CPCALC (IORDR, JKLOW(J) ,JKHIGH(J) ,JLLOW(J) ,JLHIGH(J) , 

1 JLINE ( J ) , Q , NBP , CPP) 

NBPT = NBPT + NBP 

CALL ANLIZ ( NPAT , NBP , NBPT , AREA , ANT , CPP , XC , YC , ZC) 

ENDIF 

26 CONTINUE 

DO 27 R=1 , NKSEG 

IF(KTYPE (K) .EQ. 60.0R.KTYPE (K) . EQ. 61) THEN 
IORDR=2 
NPAT=NPAT+ 1 

CALL PATCH (NORMK(K) , IORDR, KJLOW(K) , KJHIGH (K) ,KLLOW(K) , 

1 KLHIGH(K) , KLINE (K) , X,Y, Z, AREA, ANT, XC,YC,ZC) 

CALL CPCALC (IORDR, KJLOW(K) , KJHIGH (K) ,KLLOW(K) ,KLHIGH(K) , 

1 KLINE (K) , Q , NBP , CPP) 

NBPT = NBPT + NBP 

CALL ANLIZ (NPAT , NBP , NBPT , AREA, ANT , CPP , XC , YC, ZC) 

ENDIF 

27 CONTINUE 

DO 28 L=1,NLSEG 

IF(LTYPE(L) .EQ. 60.0R.LTYPE (L) .EQ. 61) THEN 
IORDR=3 
NPAT=NPAT+ 1 

CALL PATCH (NGRMTJT.) / IORDR , LJLOW ( L) , LJHIGH ( L) , LKLOW ( L) , 

1 LKHIGH(L) , LLINE (L) , X,Y, Z, AREA, ANT, XC,YC, ZC) 

CALL CPCALC ( IORDR, ULOW(L) , LJHIGH (L) ,LKLOW(L) , LKHIGH(L) , 

1 LLINE ( L) ,Q, NBP, CPP) 

NBPT = NBPT + NBP 

CALL ANLIZ (NPAT, NBP, NBPT, AREA, ANT, CPP, XC,YC,ZC) 

ENDIF 

28 CONTINUE 
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50 CONTINUE 

CALL ANLIZ ( -NPAT , NBP , NBPT , AREA , ANT , CPP , XC, YC , ZC) 

C 

•’00 FORMAT ( 1H1 , // , 5X , ' CURRENT SCALING PARAMETERS', 

1 ' FOR FORCE AND MOMENT SUMMARY' ,/,5X, 'BY PRESSURE INTEGRATION:') 
801 FORMAT (//, 5X, ' CBAR SREF SSPAN REFMX REFMY ' , 

1 ' REFMZ ' ,/ , 5X, 6F10. 3 ,/, 5X, ' RMACH ALDEG YAWDEG',/, 

2 5X, 3 F10 . 3 ) 

RETURN 

END 

SUBROUTINE PATCH (NRML, IORDR, JF, JL,KF,KL,LVL,X,Y,Z,AREA, 

* ANT, XC, YC, ZC) 

C************* ***************************** ***************************** 
cdan Routine added for force & moment calc. 7/93 

C ********************** ******************** ***************************** 

C 

include ' npami . inc ' 

DIMENSION X(JMX, KMX, LMX) , Y (JMX,KMX,LMX) ,Z(JMX,KMX,IMX) , 

1 XC(NIP) , YC(NIP) , ZC (NIP) , CPP (NIP) , AREA (NIP) ,ANT(NIP,3) 

EPSS = l.E-15 


PATCH GEOMETRY 
NP=0 

DO 100 K=KF, KL-1 
DO 100 J=JF, JL— 1 
NP=NP+1 

IF (IORDR. EQ . 1 ) THEN 
XI = X(LVL,J,K) 

Y1 » Y (LVL, J,K) 

Z1 - Z(LVL,J,K) 

X2 = X (LVL, J+l, K) 

Y2 = Y (LVL, J+l, K) 

Z2 - Z (LVL, J+l, K) 

X3 - X(LVL,J+1,K+1) 
Y3 = Y (LVL, J+1,K+1) 
Z3 = Z (LVL, J+1,K+1) 
X4 = X(LVL,J,K+1) 

Y4 = Y (LVL, J,K+1) 

Z4 = Z (LVL, J,K+1) 
ELSEIF ( IORDR . EQ . 2 ) THEN 
XI * X(J,LVL,K) 

Y1 = Y(J,LVL,K) 

Z1 = Z (J, LVL, K) 

X2 « X (J+l, LVL, K) 

Y2 = Y (J+l, LVL, K) 

Z2 *= Z (J+l, LVL, K) 

X3 =* X(J+1,LVL,K+1) 
Y3 = Y (J+l, LVL,K+1) 
Z3 » Z(J+1,LVL,K+1) 
X4 = X(J, LVL,K+1) 

Y4 - Y(J,LVL,K+1) 

Z4 = Z(J,LVL,K+1) 
ELSEIF ( IORDR . EQ . 3 ) THEN 
XI = X(J,K,LVL) 

Y1 = Y (J,K,LVL) 

Z1 = Z ( J, K, LVL) 

X2 = X (J+l , K, LVL) 
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Y2 = Y(J+1,K,LVL) 

Z2 = Z ( J+l , K, LVL) 

X3 = X(J+1,K+1,LVL) 

Y3 = Y (J+1,K+1,LVL) 

Z3 = Z (J+1,K+1,LVL) 

X4 = X(J,K+1,LVL) 

Y4 = Y (J,K+1,LVL) 

Z4 = Z ( J, K+l , LVL) 

ENDIF 

XC(NP) =. 25* (X1+X2+X3+X4) 

YC(NP) =. 25* (Y1+Y2+Y3+Y4) 

ZC(NP) =. 25* (Z1+Z2+Z3+Z4) 

XP=. 25* (X2+X3-X4— XI) 

YP=. 25* (Y2+Y3-Y4-Y1) 

ZP=. 25* (Z2+Z3-Z4-Z1) 

SMP=SQRT (XP*XP+YP*YP+ZP*ZP) 

XQ=.25*(X3+X4-X1-X2) 

YQ=. 25* (Y3+Y4-Y1-Y2) 

ZQ=. 25* (Z3+Z4-Z1— Z2) 

SMQ=SQRT (XQ*XQ+YQ*YQ+ZQ*ZQ) 

PVQX=YP*ZQ-ZP*YQ 
PVQY=ZP*XQ-XP* ZQ 
PVQZ=XP*YQ-YP*XQ 

AR=SQRT (PVQX*PVQX+PVQY*PVQY+PVQZ*PVQZ ) 

IF (AR. LT . EPSS ) AR=EPSS 
AREA (NP)=AR* 4.0 
ANT ( NP , 1 ) =PVQX/ AR 
ANT ( NP , 2 ) =PVQY/AR 
ANT (NP , 3 ) =PVQZ/AR 
IF (NRML . LT . 0 ) THEN 
ANT ( NP , 1 ) =-PVQX/AR 
ANT (NP , 2 ) =-PVQY/AR 
ANT (NP , 3 ) =-PVQZ/AR 
ENDIF 

100 CONTINUE 

RETURN 

END 

SUBROUTINE CPCALC ( IORDR, JF , JL, KF , KL, LVL , Q , NP, CPP) 

C***** *************************************** ********************** 

c 

c 

c 

include ' npami . inc ' 

COMMON/VISFAC/ ALPHA, BETA, XMACH, C2B, C2BP, VRAT, RE 
DIMENSION Q(JMX,KMX,IMX,5) ,CPP(NIP) 

GD = GAMMA *G AMI 

FSMACH = XMACH 

CPS = l./(.5*GAMMA*FSMACH**2) 

NP =0 

COMPUTE CP AT 4 'PANEL' CORNERS AND AVERAGE 

DO 5 K=KF , KL- 1 
DO 5 J=JF , JL-1 
IF(IORDR.EQ. 1)THEN 
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CALL C0RNR(LVL,J,K,CPS,GD,Q,CP1) 
CALL CORNR(LVL,J+l,K,CPS,GD,Q,CP2) 
CALL CORNR(LVL, J+l, K+l , CPS,GD, Q, CP3 ) 
CALL CORNR ( LVL, J , K+l , CPS , GD , Q , CP3 ) 
ELSEIF ( IORDR . EQ . 2 ) THEN 

CALL CORNR (J, LVL, K, CPS , GD, Q, CPI) 

CALL CORNR(J+l, LVL,K,CPS,GD, Q,CP2) 
CALL CORNR (J+l, LVL, K+l, CPS,GD,Q, CP3) 
CALL CORNR ( J , LVL , K+ 1 , CPS , GD , Q , CP4 ) 
ELSEIF ( IORDR . EQ . 3 ) THEN 

CALL CORNR (J,K, LVL, CPS , GD, Q, CPI) 
CALL CORNR(J+l,K,LVL,CPS,GD,Q,CP2) 
CALL CORNR(J+l,K+l,LVL, CPS , GD, Q, CP3) 
CALL CORNR (J , K+l , LVL, CPS , GD, Q, CP4 ) 
ENDIF 

NP = NP + 1 

CPP(NP) = 0 . 25* (CP1+CP2+CP3+CP4 ) 

5 CONTINUE 

RETURN 

END 


SUBROUTINE CORNR (J , K, L, CPS , GD, Q, CPP) 
include ' npami . inc 9 
DIMENSION Q (JMX, KMX, LMX, 5 ) 


RR 

Q1 

Q2 

Q3 

Q4 

Q5 

PP 

CPP 


l./Q(J,K,L,l) 

Q(J,K,L,1) 

Q(J,K,L,2)*RR 
Q(J,K,L,3)*RR 
Q (J, K, L, 4) *RR 
Q ( J, K, L, 5) /GAMMA 

GD* (Q5-.5* (Q2**2+Q3**2+Q4**2 ) *Q1) 
(PP-1.)*CPS 


RETURN 

END 


SUBROUTINE ANLI Z ( JP , NBP , NBPT , AREA , ANT , CPP , XC , YC , ZC ) 
include 9 npami . inc 9 
PARAMETER (MAXB=4*NIP) 

COMMON/FORCE/ CBAR , SSPAN , SREF , REFMX , REFMY , REFMZ 
COMMON/VIS FAC/ ALPHA, BETA, XMACH, C2B, C2 BP, VRAT, RE 
COMMON/SV/ WIND (3,3), SUMFX (MBC) , SUMFY (MBC) , SUMFZ (MBC) , 

1 SUMMX(MBC) ,SUMMY(MBC) , SUMMZ (MBC) 

DIMENSION CPP (NIP) , AREA (NIP) , ANT (NIP, 3) ,XC(NIP) , YC(NIP) ,ZC(NIP) 

EPS = l.E-05 
CPMIN = -10. 

CPMAX = 1.20 
JOUT = 6 
MODUP = 0 
MODWN » 0 


IF ( JP . GT . 0 ) THEN 

WRITE (JOUT , 620) JP, NBPT-NBP+1 , NBPT 
ENDIF 

IF ( JP . EQ . 1 ) THEN 
C 

C SET UP WIND AXIS SYSTEM 
C 
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OOP 


YAWR=BETA/57 . 29577951 
ALR=( ALPHA) /57. 29577951 
CB=COS (YAWR) 

CA=CB*COS (ALR) 

SA=CB*SIN (ALR) 

UXR=CA 

UYR— SIN (YAWR) 

UZR=SA 

WIND (1,1) =UXR 
WIND (2,1) =UYR 
WIND(3,1)=UZR 
WIND (1,2) =0.0 
WIND (2 , 2) =1 . 0 
WIND (3,2) =0 • 0 
WIND (1,3) =-l . 0 
WIND (2,3) =0 • 0 
WIND (3,3) =0 • 0 

IF (UZR.LT. 0. 0) WIND(1, 3) =1. 0 
US=SQRT (UXR*UXR+UYR*UYR) 

IF(US.GE. EPS) THEN 
WIND (1,2) =-UYR/US 
WIND (2,2) =UXR/US 
WIND (1,3) =-UXR*UZR/US 
WIND (2,3) =-UYR*UZR/US 
WIND (3,3) =US 
ENDIF 

ELSEIF ( JP . LT . 0 ) THEN 
WRITE (JOUT , 62 1 ) 1 , NBPT 
JP— JP 
SFX=0 . 0 
SFY-0 . 0 
SFZ=0 . 0 
SMX-0 . 0 
SMY=0 . 0 
SKZ=0 . 0 
DO 100 J=1,JP 
SFX=SFX+SUMFX ( J) /SREF 
SFY=SFY+SUMFY ( J) /SREF 
SFZ=SFZ+SUMFZ ( J) /SREF 
SMX=SMX+SUMMX ( J) /SREF 
SMY=SMY+SUMMY ( J) /SREF 
SMZ=SMZ+SUMMZ ( J) /SREF 
100 CONTINUE 

CD=SFX*WIND (1,1) +SFY*WIND (2,1) +SFZ*WIND (3,1) 

CS=S FX * WIND (1,2) +SFY*WIND (2,2) +SFZ*WIND (3,2) 

CL=S FX * WIND (1,3) +SFY*WIND (2 , 3 ) +SFZ*WIND (3,3) 

CLMX= ( SMX* WIND (1,1) +SMY*WIND (2,1) +SMZ *WIND (3,1)) /SSPAN 

CLMY= ( SMX* WIND (1,2) +SMY*WIND (2,2) +SMZ * WIND (3,2)) /CBAR 

CLMZ= (SMX* WIND (1,3) +SMY*WIND (2,3) +SMZ*WIND (3,3)) /SSPAN 

WRITE (JOUT, 629) CD, CS , CL, CIMX,CIMY , ClilZ 

SMX=SMX/ S S PAN 

SMY=SMY / CBAR 

SMZ=SMZ/SSPAN 

WRITE (JOUT , 613 ) SFX , SFY , SFZ , SMX , SMY , SMZ 
RETURN 
ENDIF 


LOADS SUMMATION 
SUMFX ( JP) =0 . 0 
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SUMFY ( JP) =0 . 0 
SUMFZ (JP)=0.0 
SUMMX ( JP) =0 . 0 
SUMMY ( JP) =0 . 0 
SUMMZ ( JP) =0 . 0 

DO 200 1=1, NBP 
XCM=XC ( I ) -REFMX 
YCM=YC ( I ) -REFMY 
ZCM=ZC ( I ) -REFMZ 
CPA=CPP ( I ) 

IF ( CPP ( I ) . LT . CPMIN ) THEN 
MODWN=MODWN+l 
CPA=CPMIN 

ELSEIF (CPP (I) . GT . CPMAX) THEN 
MODUP=MODUP+ 1 
CPA=CFMAX 
EHDIF 

DFX=-CPA* AREA ( I ) * ANT (1,1) 

DFY=-CPA* AREA ( I ) *ANT (1,2) 

DFZ=-CPA*AREA ( I ) *ANT (1,3) 

SUMFX ( JP) =SUMFX ( JP) +DFX 
SUMFY ( JP) =SUMFY ( JP) +DFY 
SUMFZ ( JP) =SUMFZ ( JP) +DFZ 
SUMMX ( JP) =SUMMX ( JP) + ( YCM*DFZ-ZCM*DFY) 

SUMMY ( JP) = SUMMY (JP) + (ZCM*DFX-XCM*DFZ) 

SUMMZ ( JP) =SUMMZ ( JP) + (XCM*DFY-YCM*DFX) 

.00 CONTINUE 

WRITE (JOUT, 699) CPMAX , MODUP , CPMIN , MODWN 

SFX-SUMFX ( JP) /SREF 

SFY-SUMFY ( JP) /SREF 

SFZ-SUMFZ ( JP) /SREF 

SMX-SUMMX ( JP) /SREF 

SMY-SUMMY ( JP) /SREF 

SMZ=SUMMZ ( JP) /SREF 

CD=SFX*WIND (1,1) +SFY *WIND (2,1) +SFZ*WIND (3,1) 

CS=SFX*WIND (1,2) +SFY*WIND (2,2) +SFZ*WIND (3,2) 

CL=SFX*WIND (1,3) +SFY*WIND (2,3) +SFZ*WIND (3,3) 

CLMX= ( SMX*WIND (1,1) +SMY*WIND (2,1) +SMZ *WIND (3,1)) /SSPAN 

CLMY= (SMX*WIND (1,2) +SMY*WIND (2,2) +SMZ*WIND(3 , 2) ) /CBAR 

CLMZ= ( SMX*WIND (1,3) +SMY*WIND (2,3) +SMZ *WIND (3,3)) /SSPAN 

WRITE (JOUT, 629) CD,CS,CL, CIMX,CIMY,CLMZ 

SMX=SMX/SS PAN 

SMY=SMY/CBAR 

SMZ=SMZ/SS PAN 

WRITE ( JOUT , 6 13 ) SFX , SFY , SFZ , SMX , SMY , SMZ 
RETURN 

605 FORMAT (// , 5X , ' INVALID WIND DIRECTION, PLEASE RE-CHECK') 

613 FORMAT ( 1H ,5X,20HIN BODY AXIS SYSTEM//9X, 3HCFX, 10X, 3HCFY, 10X, 
13HCFZ , 10X, 3HCMX, 10X, 3HCMY , 10X, 3HCMZ/1H , IX, 6G13 . 4/ , 

21H ,79 (1H*) //) 

620 FORMAT (1H0, 70 (1H*)/1H0, 

142H FORCE AND MOMENT FOR SPECIFIED PANEL SET, 12, 

2 ' ( PANELS ',15,' - ',15')',/, 

347H(NOTE: THE COEFFICIENTS HERE DO NOT INCLUDE THE, 
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412H IMAGE SET.)//) 

621 FORMAT (1H0 ,50 (1H*) /1H0 , 

123H TOTAL FORCE AND MOMENT, 

2 ' ( PANELS ',15,' - ',15')',/, 

247H(NOTE: THE COEFFICIENTS HERE DO NOT INCLUDE THE, 

312H IMAGE SET.)//) 

629 FORMAT (1H0, 

118H FORCE AND MOMENT//, 

21H , 8X, 2HCD, 11X, 2HCS , 11X, 2HCL, 10X, 3HCMX, 10X, 3HCMY, 10X, 3HCMZ/ 

31H , IX, 6G13 .4/) 

699 FORMAT (8H CPMAX= ,F6.1,18H PANELS MODIFIED= ,I5,/,8H CPMIN= ,F6.1, 
1 18H PANELS MODIFIED= ,15) 

END 
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APPENDIX B 


Program PRPAMI was developed to aid the user in pre and post-processing of data required by 
the zonal aerodynamic procedure developed under NAS3-26310. This is an interactive program written 
in FORTRAN with die specific application in mind. When practical, generality in the code logic was 
presaved, but future problems with different topological features may require software modification. 
Following is a brief description of the code features taken from die self-documentation included in die 
code. The variable IOPT is defined by the user. 


IOPT 

FEATURE 

FILES REQUIRED 

OUTPUT FILE 

1 

GENERATE NPARC RESTART 
DATA FILE 

’ .VOL’ 

REST. OLD 

2 

GENERATE VSAERO SCAN 
INPUT DATA 

INON5,REST.OLD 

SCAN.DAT 

3 

GENERATE OMNI3D GRAPHICS 
DATAFILE 

INON5, REST. OLD 

PARC.FMT 

4 

REORDER GRID BLOCKS 

’ .VOL’ 

’ .VOL2’ 

Note: Filenames in single quotes are provided by 

user . 
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onnonnnonofinnonooon to 


PROGRAM PRPAMI 


MARCH 1994 


UTILITY PROGRAM FOR DATA FILE MANIPULATION DEVEIDPED BY 
D. J. STRASH AT ANALYTICAL METHODS, INC. UNDER OPTION 2 OF 
NASA LEWIS CONTRACT NUMBER NAS3-26310. 


????QUESTIONS???? 


CALL DAN STRASH @ (206)643-9090 

OR E-MAIL dan@amiwest.com 

t ************************ ****** 


IOPT FEATURE 


FILES REQUIRED OUTPUT FILE 


1 GENERATE NPARC RESTART ' .VOL' 

DATA FILE 

2 GENERATE VSAERO SCAN IN0N5, REST. OLD 

INPUT DATA 

3 GENERATE OMNI 3D GRAPHICS IN0N5, REST. OLD 

DATA FILE 

4 REORDER GRID BLOCKS ' .VOL' 


REST. OLD 
SCAN . DAT 
PARC. FMT 
' .V0L2 ' 


(NOTE: FILENAMES IN SINGLE QUOTES ARE PROVIDED BY USER.) 


C****** ************* 


c 


CHARACTER *2 GRD 


5 CONTINUE 
WRITE ( * , * ) 
WRITE ( * , * ) 
WRITE ( * , * ) ' 
WRITE ( * , * ) 
WRITE ( * , * ) ' 
WRITE ( * , * ) ' 
WRITE ( * , * ) ' 
WRITE ( * , * ) ' 
WRITE ( * , * ) 
WRITE ( * , * ) ' 


YOU HAVE THE 

1. PREPARE 

2. PREPARE 

3. PREPARE 

4. REORDER 


PLEASE TYPE 


READ ( * , * ) IOPT 


FOLLOWING OPTIONS ' 

PARC3D RESTART FILE' 
VSAERO SCAN INPUT FILE' 
OMNI 3D PLOT FILE' 

BLOCKS (.VOL FILE)' 

NUMBER OF DESIRED OPTION. 


I 


IF (IOPT • EQ • 1) THEN 
CALL OPNRES 
CALL RVOL 
CALL PAERO 
WRITE ( * , * ) 

WRITE (*,*)' WHAT IS THE GRID TOPOLOGY?' 

WRITE (*,*)' (CO - FOR C-0 TYPE, HO - FOR H-0 TYPE)' 
WRITE (* , *) 

READ (*,80) GRD 
IF(GRD(1: 1) .EQ. 'C')THEN 
CALL COGRD 
ELSE 

CALL HOGRD 
ENDIF 

ELSEIF ( IOPT . EQ . 2 ) THEN 
CALL OPNSCN 
CALL READRES 
CALL SCNGEN 
ELSEIF ( IOPT . EQ . 3 ) THEN 
CALL OPNPLT 
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CALL READRES 
CALL PLTGEN 
ELSEIF ( IOPT . EQ . 4 ) THEN 
CALL OPNRDR 
CALL REORDR 
ELSEIF ( IOPT . LT . 0 ) THEN 
STOP 
ELSE 

WRITE (*,*)' UNSUPPORTED OPTION, TRY AGAIN' 

GO TO 5 
ENDIF 
STOP 

80 FORMAT (A2) 

END 

SUBROUTINE OPNRES 

CHARACTER* 80 FILUN, FI LOT, FILIN 
DATA FILOT/ 'REST. OLD'/ 

WRITE ( * f * ) 9 INPUT NAME OF GRIDGEN VOLUME GRID FILE.' 

READ(*, 80) FILIN 

OPEN (UNIT 3 * 8 , FILE-FILIN, FORM* 'UNFORMATTED' , STATUS 3 *' OLD' ,err=9 001) 
OPEN (UNIT=9 , FILE=FILOT , FORM= 'UNFORMATTED' , STATUS 3 * ' UNKNOWN ' , 

1 err=9002) 

REWIND 8 
REWIND 9 
RETURN 

9001 WRITE (*,*)' UNABLE TO OPEN ', FILIN 
STOP 

;Q2 WRITE (*, *) ' UNABLE TO OPEN REST. OLD' 

STOP 

80 FORMAT (A80) 

END 

SUBROUTINE OPNSCN 

CHARACTER* 80 FILUN, FILOT, FILIN 

DATA FILIN/ 'INON5'/, FILUN/ 'REST . OLD'/ , FILOT/ ' SCAN . DAT ' / 

OPEN (UNIT=7 , FILE=FILIN , FORM— ' FORMATTED ' , STATUS** ' OLD ' , err=9 001 ) 
OPEN (UNIT=8 , FILE 3 * FILUN , FORM** ' UNFORMATTED ' , STATUS 3 * ' OLD ' , err=9 002 ) 
OPEN (UNIT=9 , FILE=*FILOT , FORM* 3 ' FORMATTED ' , STATUS** 'UNKNOWN ' , 

1 err*=9003) 

REWIND 7 
REWIND 8 
REWIND 9 
RETURN 

9001 WRITE ( * f * ) ' UNABLE TO OPEN INON5' 

STOP 

9002 WRITE (*, *) ' UNABLE TO OPEN REST. OLD' 

STOP 

9003 WRITE (* , *) ' UNABLE TO OPEN SCAN , DAT' 

STOP 

END 

SUBROUTINE OPNPLT 

CHARACTER* 80 FILUN, FILOT, FILIN 

DATA FILIN/ 'INON5'/, FILUN/ 'REST. OLD'/ , FILOT/ 'PARC. FMT'/ 
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OPEN(UNIT=7,FILE=FILIN,FORM='FORMATTED' ,STATUS='OLD' ,err-9001) 

OPEN (UNIT=8 , FILE=FILUN f FORM= ' UNFORMATTED ' , STATUS='OLD' ,err=9002) 
OPEN (UNIT=9 , FILE=FILOT , FORM= ' FORMATTED ' , STATOS= ' UNKNOWN ' , err=9 003) 

REWIND 7 
REWIND 8 
REWIND 9 
RETURN 

9001 WRITE (*,*)' UNABLE TO OPEN INON5' 

STOP 

9002 WRITE (* , *) ' UNABLE TO OPEN REST. OLD' 

STOP 

9003 WRITE ( * , * ) ' UNABLE TO OPEN PARC . FMT ' 

STOP 

END 

SUBROUTINE OPNRDR 

CHARACTER* 8 0 FILUN, FI LOT, FILIN 

DATA FI LOT/' REST. OLD'/ 

WRITE (*,*)' INPUT NAME OF GRIDGEN VOLUME GRID FILE. ' 

READ (*,80) FILIN 

OPEN (UNIT=8 , FILE=FILIN , FORM= ' UNFORMATTED ' ,STATUS='OLD' ,err=9001) 
WRITE (*,*)' INPUT NAME OF MODIFIED GRIDGEN VOLUME GRID FILE.' 
READ(* ,80) FILOT 

OPEN (UNIT=9 , FILE=FILOT , FORM= 'UNFORMATTED ' , STATUS= 'UNKNOWN ' , 

1 err-9002) 

REWIND 8 
REWIND 9 
RETURN 

9001 WRITE ( * , * ) ' UNABLE TO OPEN ', FILIN 
STOP 

9002 WRITE ( * , * ) ' UNABLE TO OPEN FILOT 
STOP 

80 FORMAT (A80) 

END 

subroutine rvol 
c Read the VOLGRID file. 

PARAMETER (NTB=3, JMX=230, KMX=32, LMX=50) 

common /grid/ x(jmx,kmx,lmx,NTB) ,y(jmx,kmx,lmx,NTB) , 

& z (jmx,kmx,lmx,NTB) , jm(NTB) , km(NTB) ,lm(NTB) ,NBIOCKS,NBB,SYM 


io3di = 8 


c Read the volgrid file in gridgen3d format 

WITHOUT the ghost cells. 


read ( io3di , end=9010 , err=902 0 ) nblocks 
IF ( NBLOCKS . GT . NTB ) THEN 
writ»{*. Toil) nblocks 
ELSE 

write (*,1030) nblocks 
ENDIF 

do 100 mb = 1, nblocks 
read(io3di) jm(mb), km(mb) , lm(mb) 


51 



write(*, 1040) mb, jm(mb) , Jan(mb) , lm(mb) 


& 

& 

& 

& 

& 

& 

6 

& 

& 

& 

& 


read(io3di,end=9010,err=9020) ( ( ( x(i, j ,k,mb) , 

i = 1, jm(mb) ) , 

j = lr Jan (Mb) ) , 

k - 1, lm(mb) ) , 

( ( ( y(i, j,k,mb) , 

i = 1, jm(mb) ) , 
j = 1, km (mb) ) , 
k - 1 , lm(mb) ) , 

( ( ( z(i, j,k,mb) f 

i = 1, jm(mb) ) , 
j = 1, km(rnb) ) , 
k = 1, lm(mb) ) 


100 continue 


close (unit=io3di) 
return 

9010 write (*, 1010) 
stop 

9020 write (*, 1020) 
stop 

1010 format (/, 5X, ' ERROR: END OF FILE READING V0LGRID FILE', 

& /) 

"20 format (/ , 5X , ' ERROR : ERROR READING VOLGRID FILE', 

& /) 

1030 format (/, 5X, ' NUMBER OF BLOCKS IN VOLGRID FILE... '^3, 

* /) 

1031 format(/,5X, ' ERROR READING VOLGRID FILE...',/, 

& 5X, ' INCREASE BLOCK DIMENSIONS (NTB) TO ',i3, 

5 /) 

1040 format (/,5X, ' NO. OF PTS IN BLOCK ',i3,' OF VOLGRID FILE', 

6 /,5X, ' VOLGRID - ' , i3, ' X ',i3,' X ',i3, 

& /) 

end 

subroutine reordr 
c Read the VOLGRID file. 

PARAMETER (NTB=3, JMX=230, KMX=32, UtX=50) 

ccanmon /grid/ x(jmx,Janx,lmx,NTB) ,y(jmx,Janx,lmx,NTB) , 

& z ( jmx,Janx, lmx,NTB) , jm(NTB) , km (NTB) ,1m (NTB) , NBLOCKS , NBB , SYM 
DIMENSION IB (NTB) 


io3di = 8 
io3do = 9 


c Read the volgrid file in gridgen3d format 

c WITHOUT the ghost cells. 


read ( io3di , end=9010 , err=9020) nblocks 
write (io3do) nblocks 
write (*,1030) nblocks 

do 100 mb = 1, nblocks 
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& 

& 

& 

& 

& 

& 

& 

& 

& 

& 

& 


read(io3di) jm(mb) , lon(inb) , lm(mb) 
write (*,1040) mb, jm(mb) ,k 2 n(mb) 

read ( io3di , end=9010 , err=902 0 ) ( ( ( x(i, j ,k,inb) i 

i=l, jm(mo) ) , 
j = 1, km (mb) ) , 
k = 1, lm(mb) ) , 

( ( ( y(i,j,k,mb), 

i = 1, jm(mb) ) , 
j = 1, km(mb) ), 
k = 1, lm(mb) ) , 

( ( ( z(i,j,k,mb), 

i = 1, jm(mb) ) , 
j = 1, km (mb) ) , 
k = 1, lm(mb) ) 


100 continue 


££ite(*!*)' WOULD YOU LIKE TO SCALE AND TRANSLATE GRID?' 
write(*,*)' (1 = YES, 0 = NO)' 

read ( * , * ) iopt 
if ( iopt . ne . 0 ) then 

write ( * , * ) ' INPUT REFERENCE LENGTH , XT , YT , ZT ' 
write ( * , * ) 

read ( * , * ) refl,xt,yt,zt 
if (abs(refl) .It. l.e-07) then 
write ( * , * ) 

write (*,*) ' TRY AGAIN...' 
go to 5 
endif 

do 150 mm=l, nblocks 
do 155 j=l, jm(mm) 
do 155 k=l,km(mm) 
do 155 1=1, lm(mm) 

x(j,k,l,mm)=x(j,k,l,mm)/refl + xt 
y(j,k,l,mm)=y(j,k,l,mm)/refl + yt 
z(j,k,l,mm)=z(j,k,l,mm)/refl + zt 

155 continue 
150 continue 
endif 

write i *) 

write(*,*)' PLEASE TYPE IN NEW BIOCK ORDER, ' 

write (*,*)' (e-9* 1 » 3 , 2 ) . ' 

write (*, *) 

read(*,*) (ib(n) ,n=l, nblocks) 
do 200 mm = 1, nblocks 


& 

& 

& 

& 


mb = ib(mm) . 

write (io3do) jm(mb) , km (mb) , Iststf) 
write (*, 1040) mb, jm(mb) , tan (mb) ,lm(mb) 


write (io3do) 


( ( ( x(i, j ,k,mb) , 

i = 1 , jm(mb) ) , 
j = 1, km(mb) ) , 
k = 1, lm(mb) ) , 

( ( ( Y(i/ j / 
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& 

& 

& 

& 

& 

& 

& 

200 continue 


close (unit=io3di) 
close (unit=io3do) 
return 

9010 write (*, 1010) 
stop 

9020 write (*, 1020) 
stop 

1010 format (/ , 5X, 'ERROR: END OF FILE READING VOLGRID FILE*, 

* /) 

1020 format (/,5X, 'ERROR: ERROR READING VOLGRID FILE', 

& /) 

1030 format (/ , 5X, ' NUMBER OF BLOCKS IN VOLGRID FILE ',i3, 

& /) 

1040 format (/ , 5X , ' NO. OF PTS IN BLOCK ',i3,' OF VOIGRID FILE', 

& / , 5X, ' VOLGRID - ' ,i3, ' X ',i3,' X ',i3, 

t /) 

end 

SUBROUTINE PAERO 

PARAMETER (NTB=3, JMX=230, KMX=32, IMX*50) 

common /grid/ x(jmx,kmx,lmx,NTB) ,y ( jmx,kmx, lmx,NTB) , 

& z(jmx,kmx,lmx,NTB) , jm(NTB) , km(NTB) ,lm(NTB) , NBLOCKS , NBB, SYM 
COMMON/AERO/ R( jmx,kmx, lmx) ,RU( jinx, kmx, lmx) ,RV( jmx,kmx,lmx) , 

1 RW( jmx,kmx, lmx) ,E( jmx,kmx, lmx) , ALPHA , BETA , FMACH , CBAR 
CHARACTER *1 GOM 
LOGICAL SYM 

DATA G/1.4/, GM1/0.4/, SYM/. FALSE./ 

WRITE(*,*) 

WRITE (*,*)' INPUT ALPHA, BETA, MACH NUMBER AND CBAR...' 

WRITE (*,*)' (CBAR USED TO SCALE GEOMETRY)' 

WRITE ( * * ) 

READ (*,*) ALDEG , YAWDEG , FMACH , CBAR 
IF ( ABS (YAWDEG) . GE . 0 . 1) THEN 
WRITE ( * / * ) 

WRITE (*!*)' SHOULD THE GEOMETRY BE MIRRORED ACROSS X-Z PLANE?' 
WRITE (*, *) ' (T-TRUE, F- FALSE) ' 

WRITER,*) 

READ (*,80) GOM 

IF(GOM(l:l) .EQ. ' T ' ) SYM= . TRUE . 

ENDIF 

IF ( FMACH . LT . 0 . 1 ) FMACH=0 . 1 
NBB-NBLOCKS 
IF(SYM) NBB=»2*NBLOCKS 
IF (NBLOCKS . GT . NTB) THEN 


i * 1, jm(mb) ) , 
j * 1, Xm(mb) ) , 
k * 1, lm(mb) ) , 

( ( ( z(i, j,k,mb) , 

i = 1, jm(mb) ) , 
j = 1, km(mb) ) , 
k = 1, lm(mb) ) 
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write (*,1031) nblocks 
ENDIF 


• FORM THE ARRAYS OF NON-DIMENSIONAL CONSERVATION VARIABLES 
C CONSISTENT WITH THE FREE -STREAM MACH NUMBER 

Y AWR= YAWDEG/ 5 7 . 29577951 
ALR = ALDEG/57. 29577951 
CB=COS(YAWR) 

CA=CB*COS ( ALR) 

SA=CB*SIN (ALR) 

FACT= ( 1 . + . 2*FMACH**2 ) 

PBAR=FACT** (-3 . 5) /G 
RHOB=FACT** (-2 .5) 

RHOU=RHOB*FMACH*SQRT ( 1 . /FACT) 

EBAR=PBAR/GM1+ . 5* (RHOU**2 ) /RHOB 
RHOU =RHOU * CA 
RHOV=-RHOU*SIN (YAWR) 

RHOW=RHOU * S A 

DO 4 L=1 , IMX 
DO 4 K=1,KMX 
DO 4 J=l, JMX 
R (J,K,L)=RHOB 
RU(J,K,L) =RHOU 
RV(J,K,L)=RHOV 
RW(J,K,L)=RHOW . 

E (J,K,L)=EBAR 
/ CONTINUE 

return 

80 FORMAT (Al) 

81 FORMAT (F10.0) . , 

1031 format (/,5X, ' ERROR READING VOLGRID FILE... ,/, 

& 5X, ' INCREASE BLOCK DIMENSIONS (NTB) TO ',i3, 

& /) 

END 


SUBROUTINE HOGRD 

C***** ************************ 

C ASSUMES THE FOLLOWING: 

C 1. 3-BLOCK H-0 TOPOLOGY ^ . , , 

C 2. BLOCKS AND INDICES ARE IN PROPER ORDER (upper-tip- lower 

C* ********************************************* ******** ******** 

c 

PARAMETER (NTB=3, JMX=230, KMX=32, LMX=50) 

common /grid/ x (jmx, kmx, lmx, NTB) ,y (jmx,kmx,lmx,NTB) , 

& z( jmx, kmx, lmx, NTB) , jm(NTB) , km (NTB) , lm (NTB) , NBLOCKS, NBB , SYM 
COMMON/AERO/ R( jmx, kmx, lmx) ,RU( jmx, kmx, lmx) ,RV( 3 mx,kmx,lmx) , 
1 RW( jmx, kmx, lmx) ,E(jmx, kmx, lmx) , ALPHA, BETA, FMACH,CBAR 
DIMENSION XX (JMX, KMX, LMX) ,YY (JMX, KMX, IMX) , ZZ (JMX, KMX, IMX) 
LOGICAL SYM 
DATA G/l. 4/ , GM1/0.4/ 


jto = 9 

WRITE RESTART TAPE 
NC1=0 

WRITE (jto) NCI, G 

WRITE (*,*) 'JMAX KMAX LMAX' 
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DO 2 MB= 1 , NBLOCKS 
C BLOCK #1: STARBOARD UPPER 
IF (MB . EQ . 1 ) THEN 

JMAX-JM(MB) 

KMAX— KM (MB) 

LMAX— IM (MB) +1 
WRITE(*,*) JMAX, KMAX, IMAX 
DO 5 J-l, JMAX 
DO 5 K=1 , KMAX 
DO 5 L— 1 , LMAX 
IF(L.NE. 1) THEN 

XX(J,K,L)=X(J,K,L-1,MB)/CBAR 

YY(J,K,L)-Y(J,K,L-1,MB)/CBAR 

ZZ(J,K,L)=Z(J,K,L-1,MB)/CBAR 

ELSE 

XX (J,K,L)=X(J, KM (NBLOCKS) -K+l, L+l, NBLOCKS) /CBAR 
YY ( J , K, L) =Y ( J , KM (NBLOCKS ) -K+l , L+l , NBLOCKS) /CBAR 
ZZ ( J , K, L) =Z ( J , KM (NBLOCKS ) -K+l , L+l , NBLOCKS) /CBAR 
ENDIF 
5 CONTINUE 

CALL WRDAT(JTO, JMAX, KMAX, LMAX, XX, YY, ZZ ,R,RU,RV,RW,E) 

ELSEIF (MB . EQ . NBLOCKS ) THEN 

C BLOCK 3: STARBOARD LOWER 
JMAX-JM(MB) 

KMAX— KM (MB) 

IMAX— IM (MB) 

WRITE (*,*) JMAX , KMAX , IMAX 
DO 15 J=l, JMAX 
DO 15 K=1 , KMAX 
DO 15 L*1,IMAX 
XX (J r K, L) =X( J ,K, L, MB) /CBAR 
YY ( J , K , L ) =Y ( J , K , L , MB ) / CBAR 
ZZ (J,K,L) =Z ( J , K, L, MB) /CBAR 
15 CONTINUE 

CALL WRDAT(JTO, JMAX, KMAX, LMAX, XX, YY,ZZ,R,RU,RV,RW,E) 
ELSE 

C BLOCK #2: STARBOARD TIP 
JMAX=JM(MB) 

KMAX=KM(MB)+2 
IMAX— LK(MB) 

WRITE (*,*) JMAX , KMAX , LMAX 
DO 10 J=l, JMAX 
DO 10 K=1 , KMAX 
DO 10 L»l, IMAX 
IF(K.EQ. 1)THEN 

XX ( J , K, L) =X ( J , KM(MB-l) -1 , L, MB“-1) /C2AR 
YY (J, K, L) *Y (J,KM(MB-1) -1, L,MB-1)/CBAR 
ZZ ( J , K, L) =Z ( J , KM (MB-1 ) -1 , L, MB-1) /CBAR 
ELSEIF ( K . EQ . KMAX ) THEN 

XX(J,K,L)=X(J,2,L,MB+1)/CBAR 
YY ( J , K, L) =Y ( J , 2 , L, MB+1 ) /CBAR 
ZZ ( J , K, L) =Z ( J , 2 , L, MB+1) /CBAR 
ELSE 
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XX(J,K,L)=X(J,K-1,L,MB)/CBAR 

YY(J,K,L)=Y(J,K-1,L,MB)/CBAR 

ZZ ( J , K, L) =Z ( J , K-l , L,MB) /CBAR 
ENDIF 

10 CONTINUE 

CALL WRDAT(JTO, JMAX , KMAX , LMAX , XX , YY , ZZ ,R,RU,RV,RW,E) 
ENDIF 

2 CONTINUE 

IF ( . NOT . S YM) GO TO 99 
DO 3 MB=NBLOCKS, 1,-1 
IF (MB . EQ . NBLOCKS ) THEN 


C 


BLOCK 4: PORT LOWER 
JMAX=JM (MB) 

KMAX=KM (MB) +1 
LMAX=LM(MB) 

WRITE (*,*)JKAX, KMAX, LMAX 

DO 20 J=l, JMAX 
DO 20 K=l, KMAX 


DO 20 L=l, LMAX 
IFfK.GT. 1)THEN 

XX(J, K,L) = X( J , KM(MB) -K+2 , L,MB) /CBAR 
YY ( J , K, L) =-Y ( J , KM (MB) -K+2 , L , MB) /CBAR 
ZZ ( J, K, L) = Z (J,KM(MB) -K+2 ,L, MB) /CBAR 


ELSE 

XX(J,K,L) = X( J,KM(MB) — 1, L,MB) /CBAR 
YY ( J , K, L) = Y (J,KM(MB) — 1 , L,MB) /CBAR 
ZZ ( J, K, L) = Z ( J ,KM(MB) -1 , L,MB) /CBAR 
ENDIF 

20 CONTINUE 


CALL WRDAT(JTO, JMAX, KMAX, LMAX, XX, YY, ZZ ,R,RU,RV,RW,E) 


ELSEIF (MB . EQ . 1) THEN 


C 


BLOCK #6: PORT UPPER 
JMAX=JM(MB) 

KMAX=KM(MB)+1 
1MAX=LM(MB) +1 
WRITE (*,*) JMAX , KMAX , IMAX 


DO 30 J=l, JMAX 
DO 30 K=1 , KMAX 
DO 30 1^1, LMAX 
IF ( K . NE . KMAX ) THEN 
IF(L.NE.1)THEN 
XX ( J , K, L) = X(J, 
YY (J , K, L) =-Y ( J , 
ZZ(J,K,L)= Z(J, 
ELSE 

XX ( J , K, L) = X(J, 
YY (J,K,L)=-Y (J , 
ZZ (J,K,L)= Z(J, 


KM (MB) -K+l , L-l , MB) /CBAR 
KM (MB) -K+l , L-l ,MB) /CBAR 
KM (MB) -K+l , L-l , MB) /CBAR 

K,L+1, NBLOCKS) /CBAR 
K, L+l , NBLOCKS) /CBAR 
K,L+1, NBLOCKS) /CBAR 


ENDIF 

ELSE 

IF (L.NE . 1) THEN 

XX ( J, K, L) = X(J, 2 , L-l, MB) /CBAR 
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YY(J,K,L)= Y (J, 2 , L-l ,MB) /CBAR 
ZZ (J, K,L) = Z (J, 2, L-l, MB) /CBAR 
ELSE 

XX ( J , K, L) = X ( J , KM ( NBLOCKS ) -1, L+l,NBLOCKS)/CBAR 
YY(J,K,L)= Y (J, KM (NBLOCKS )-l,L+l, NBLOCKS) /CBAR 
ZZ (J, K, L) = Z (J, KM(NBLOCKS) -1, L+l, NBLOCKS) /CBAR 
ENDIF 
ENDIF 

30 CONTINUE 

CALL WRDAT ( JTO , JMAX , KMAX , IMAX , XX , YY , ZZ , R, RU , RV , RW, E) 
ELSE 

C BLOCK #5: PORT TIP 
JMAX=JM(MB) 

KMAX— KM (MB) +2 
LMAX=LM(MB) 

WRITE (*,*) JMAX, KMAX, LMAX 
DO 25 J=l, JMAX 
DO 25 K=1 , KMAX 
DO 25 L=1 , LMAX 
IF (K. EQ . KMAX) THEN 

XX ( J , K, L) = X( J , KM( 1) -1, L, 1) /CBAR 
YY (J ,K, L) =-Y ( J , KM( 1) -1, L, 1) /CBAR 
ZZ(J,K,L)= Z ( J , KM( 1) -1 , L, 1) /CBAR 
ELSEIF ( K . EQ , 1 ) THEN 

XX ( J , K, L) = X (J, 2 ,L, NBLOCKS) /CBAR 
YY (J , K , L) =-Y (J , 2 , L, NBLOCKS ) /CBAR 
ZZ(J,K,L)= Z (J, 2, L, NBLOCKS) /CBAR 
ELSE 

XX ( J , K, L) = X(J ,KM(MB) — K+2 , L,MB)/CBAR 
YY ( J, K, L) =-Y ( J, KM(MB) -K+2 , L, MB) /CBAR 
ZZ ( J, K,L) = Z ( J, KM (MB) -K+2, L, MB) /CBAR 
ENDIF 

25 CONTINUE 

CALL WRDAT (JTO, JMAX, KMAX, LMAX, XX, YY,ZZ,R,RU,RV,RW,E) 

ENDIF 
3 CONTINUE 

99 RETURN 

9001 write(*, 1050) 
stop 

1050 format (/,5X, 'ERROR OPENING OUTPUT FILE 

& /) 


END 

SUBROUTINE COGRD 

****************** 

C ASSUMES THE FOLLOWING: 

1. 3 -BLOCK C-0 TOPOLOGY 
C 2. BLOCKS AND INDICES ARE IN PROPER ORDER 
c*** ********* ****** 

c 

PARAMETER (NTB=3 , JMX=230, KMX=32, LMX=50) 
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common /grid/ x(jmx,kmx,lmx,NTB) ,y(jmx,kmx,lrax,NTB) , 

& Z ( jmx, kmx, lmx, NTB) , jm(NTB) , km (NTB) , lm(NTB) , NBIOCKS ,NBB, SYM 
COMMON/AERO/ R( jmx, kmx, lmx) ,RU (jmx, kmx, lmx) ,RV(jmx,kmx,lmx) , 
1 RW (jmx, kmx, lmx) ,E (jmx, kmx, lmx) , ALPHA , BETA , FMACH , CBAR 
DIMENSION XX (JMX, KMX, LMX) , YY (JMX, KMX, LMX) , ZZ (JMX, KMX, LMX) 
LOGICAL SYM 
DATA G/1.4/, GM1/0.4/ 

jto = 9 

C WRITE RESTART TAPE 
NC1=0 

WRITE (jto) NCI, G 

WRITE (*,*) 'JMAX KMAX LMAX ' 

DO 2 MB=1 , NBLOCKS 
C BLOCK #1: STARBOARD UPPER 
IF (MB . EQ . 1 ) THEN 

JMAX=JM ( MB ) + 1 
KMAX=KM(MB) 

LMAX=LM (MB) +1 
WRITE (*,*) JMAX, KMAX, LMAX 
DO 5 J=2, JMAX 
DO 5 K=1,KMAX 
DO 5 L=1 , LMAX 
IF ( L . NE . 1 ) THEN 

XX(J,K,L)=X(J-1,K,L-1,MB)/CBAR 
YY ( J , K, L) =Y ( J.-l , K, L-l , MB) /CBAR 
ZZ ( J,K, L) =Z ( J-l ,K, L-1,MB) /CBAR 
ELSE 

XX ( J , K , L) =X (J-l , KM (NBLOCKS ) -K+l , L+l , NBLOCKS) /CBAR 
YY ( J , K, L) =Y (J-l - KM (NBLOCKS ) -K+l , L+l , NBLOCKS) /CBAR 
ZZ ( J , K, L) =Z (J-l , KM (NBLOCKS ) -K+l , L+l , NBLOCKS) /CBAR 
ENDIF 
5 CONTINUE 

DO 8 K=1 , KMAX 
DO 8 L=l, LMAX 
IF(L.NE.1)THEN 

XX (1, K, L) =X (2 , KM (NBLOCKS) -K+l, L-l, NBLOCKS) /CBAR 
YY ( 1 , K , L) =Y (2 , KM (NBLOCKS ) -K+l , L-l , NBIOCKS ) /CBAR 
ZZ (1,K,L)=Z (2 , KM (NBLOCKS) -K+l, L-l, NBIOCKS) /CBAR 
ELSE 

XX ( 1 , K , L) =X ( 3 , KM (NBLOCKS ) -K+l , L , NBLOCKS ) /CBAR 
YY ( 1 , K, L) =Y ( 3 , KM (NBLOCKS ) -K+l , L, NBIOCKS ) /CBAR 
ZZ ( 1 , K, L) =Z ( 3 , KM (NBLOCKS ) -K+l , L, NBIOCKS) /CBAR 
ENDIF 
8 CONTINUE 

CALL WRDAT ( JTO, JMAX, KMAX, LMAX, XX, YY, ZZ ,R,RU,RV,RW,E) 

ELSEIF (MB . EQ . NBIOCKS ) THEN 

C BLOCK 3: STARBOARD LOWER 
JMAX=JM (MB) 

KMAX=KM (MB) 

IMAX=*LM(MB) 

WRITE (*,*) JMAX, KMAX, LMAX 
DO 15 J=l, JMAX 
DO 15 K=1 , KMAX 
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DO 15 L=l, LMAX 
XX(J,K,L)=X(J,K,L,MB)/CBAR 
YY(J,K,L)=Y(J,K,L,MB)/CBAR 
ZZ(J , K,L) =Z (J,K,L,MB)/CBAR 
15 CONTINUE 

CALL WRDAT(JTO, JMAX, KMAX, IMAX, XX, YY,ZZ,R,RU,RV,RW,E) 
ELSE 

C BLOCK #2: STARBOARD TIP 
JMAX=JM(MB) 

KMAX=KM(MB) +2 
LMAX=LM(MB) 

WRITE (*,*) JMAX, KMAX, LMAX 
DO 10 J*=l, JMAX 
DO 10 K=1,KMAX 
DO 10 L=l, LMAX 
IF(K.EQ. 1) THEN 

XX ( J , K, L) =X ( J , KM (MB-1) -1 , L, MB-1) /CBAR 
YY(J ,K,L) =Y ( J , KM (MB-1) -1, L, MB-1) /CBAR 
ZZ (J, K, L) =Z (J , KM (MB-1 ) -1 , L, MB-1 ) /CBAR 
ELSEIF ( K . EQ . KMAX ) THEN 

XX(J,K,L)=X(J,2,L,MB+1)/CBAR 
YY(J,K,L)=Y(J,2,L,MB+1)/CBAR 
ZZ ( J,K, L) =Z ( J, 2 , L,MB+1) /CBAR 
ELSE 

XX(J,K,L)~X(J,K-1,L,MB)/CBAR 
YY ( J, K, L) =Y ( J , K-l , L,MB) /CBAR 
ZZ ( J , K, L) =Z ( J , K-l , L, MB) /CBAR 
ENDIF 

10 CONTINUE 

CALL WRDAT(JTO, JMAX, KMAX, LMAX, XX, YY, ZZ,R,RU r RV,RW,E) 

ENDIF 
2 CONTINUE 

IF( .NOT.SYM) GO TO 99 

DO 3 MB=NBLOCKS ,1,-1 

IF (MB . EQ . NB LOCKS ) THEN 

C BLOCK 4: PORT LOWER 
JMAX=JM (MB) 

KMAX=KM (MB) +1 
LMAX=LM (MB) 

WRITE ( * , * ) JMAX , KMAX , LMAX 
DO 20 J=l, JMAX 
DO 20 K=1 , KMAX 
DO 20 L=1 , LMAX 
IF(K.GT.1)THEN 

XX ( J, K, L) = X ( J, KM (ME) -K+2 , L, MB) /CBAR 
YY ( J , K, L) =- Y ( J , KM (MB) -K+2 , L, MB) /CBAR 
ZZ (J , K, L) = Z (J, KM (MB) -K+2 ,L, MB) /CBAR 
ELSE 

XX(J,K, L) = X(J,KM(MB) -1,L,MB) /CBAR 
YY ( J, K, L) = Y ( J , KM(MB) -1, L,MB) /CBAR 
ZZ ( J , K, L) = Z ( J , KM(MB) -1, L,MB) /CBAR 
ENDIF 
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20 CONTINUE 


CALL WRDAT ( JTO , JMAX , KMAX , LMAX , XX, YY , ZZ , R, RU , RV, RW, E) 
ELSEIF (MB . EQ . 1 ) THEN 

C BLOCK #6: PORT UPPER 
JMAX=JM (MB) +1 
KMAX=KM(MB)+1 
IMAX=LM(MB)+1 
WRITE (*,*) JMAX , KMAX , LMAX 
DO 30 J=2 , JMAX 
DO 30 K=1 , KMAX 
DO 30 L^l , LMAX 
IF ( K . LT . KMAX ) THEN 
IF(L.NE.1)THEN 

XX(J,K,L) =X( J-l, KM (MB) -K+l, L-l, MB) /CBAR 
YY ( J , K, L) =-Y ( J-l , KM (MB) -K+l , L-l , MB) /CBAR 
ZZ (J, K, L) =Z (J-l, KM (MB) -K+l, L-l, MB) /CBAR 
EESE 

XX ( J, K,L)=X( J-l, K, L+l, NBLOCKS) /CBAR 
YY ( J , K , L) =-Y ( J- 1 , K , L+ 1 , NBLOCKS ) /CBAR 
ZZ ( J , K, L) =Z ( J-l , K, L+l , NBLOCKS) /CBAR 
ENDIF 
ELSE 

IF(L.NE. 1)THEN 

XX ( J , K, L) =X ( J-l , 2 , L-l , MB) /CBAR 
YY (J, K, L) =Y ( J;-l , 2 , L-l , MB) /CBAR 
ZZ ( J , K, L) =Z ( J-l , 2 , L-l , MB) /CBAR 
ELSE 

XX(J,K,L) =X ( J-l , KM (NBLOCKS) —1, L+l , NBLOCKS) /CBAR 
YY ( J , K, L) =Y ( J-l , KM (NBLOCKS ) -1 , L+l , NBLOCKS) /CBAR 
ZZ ( J , K, L) =Z ( J-l , KM (NBLOCKS ) -1 , L+l , NBLOCKS ) / CBAR 
ENDIF 
ENDIF 

30 CONTINUE 

DO 32 K=1 , KMAX 
DO 32 L=1 , LMAX 
IF ( K - LT . KMAX ) THEN 
IF(L.NE.1)THEN 

XX ( 1 , K, L) =X ( 2 , K , L-l , NBLOCKS) /CBAR 
YY ( 1 , K, L) =-Y ( 2 , K, L-l , NBLOCKS) /CBAR 
ZZ ( 1 , K , L) =Z ( 2 , K , L-l , NBIOCKS ) / CBAR 
ELS£ 

XX ( 1 , K, L) =X ( 3 , K, L, NBLOCKS ) /CBAR 
YY ( 1 , K, L) =- Y ( 3, K,L, NBIOCKS) /CBAR 
ZZ (1,K, L) =Z ( 3, K,L, NBLOCKS) /CBAR 
ENDIF 
ELSE 

IF(L.NE.1)THEN 

XX ( 1 , K, L) =X ( 2 , KM (NBLOCKS) -1 , L-l, NBLOCKS) /CBAR 
' YY ( 1, K, L) =Y (2 , KM (NBLOCKS) —1, L—l, NBIOCKS) /CBAR 

ZZ ( 1 , K , L) =Z ( 2 , KM (NBIOCKS ) -1 , L-l , NBIOCKS ) /CBAR 

ELSE 

XX ( 1 , K , L) =X ( 3 , KM ( NBLOCKS ) -1 , L , NBIOCKS ) /CBAR 
YY ( 1 , K , L) =*Y ( 3 , KM ( NBLOCKS )- 1 , L , NBIOCKS ) /CBAR 
Z Z ( 1 , K , L ) =Z ( 3 , KM ( NBLOCKS )- 1 , L , NBIOCKS ) /CBAR 
ENDIF 
ENDIF 
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32 CONTINUE 


CALL WRDAT (JTO, JMAX, KMAX, LMAX, XX, YY, ZZ , R,RU,RV,RW,E) 
ELSE 

C BLOCK #5: PORT TIP 
JMAX=JM (MB) 

KMAX=KM (MB) +2 
LMAX=IM (MB) 

WRITE ( * , * ) JMAX , KMAX , LMAX 
DO 25 J=l, JMAX 
DO 25 K=1 , KMAX 
DO 25 L*=1,LMAX 
IF ( K . EQ . KMAX ) THEN 

XX ( J , K, L) = X( J, KM ( 1) ~1 , L, 1) /CBAR 
YY (J,K,L) =-Y ( J, KM (1) -1, L, 1)/CBAR 
ZZ (J,K,L) = Z ( J, KM (1) -1 , L, 1) /CBAR 
ELSEIF ( K . EQ . 1 ) THEN 

XX ( J , K, L) = X(J, 2,L,NBLOCKS)/CBAR 
YY(J,K,L)=-Y(J,2,L,NBLOCKS)/CBAR 
ZZ(J,K,L)= Z (J, 2, L, NBLOCKS) /CBAR 
ELSE 

XX (J , K, L) = X ( J , KM (MB) -K+2 , L, MB) / CBAR 
YY ( J , K, L) =-Y ( J , KM (MB) -K+2 , L, MB) /CBAR 
ZZ (J, K, L) = Z (J r KM(MB) -K+2 , L, MB) /CBAR 
ENDIF 

25 CONTINUE 

CALL WRDAT (JTO, JMAX, KMAX, LMAX, XX, YY,ZZ,R,RU,RV,RW,E) 

ENDIF 
3 CONTINUE 

99 RETURN 

9001 write (*, 1050) 

Stop 

1050 format(/,5X, 'ERROR OPENING OUTPUT FILE ', 

6 /) 


END 

SUBROUTINE WRDAT(JTO, JMAX, KMAX, IMAX, XX, YY,ZZ,R,RU,RV,RW,E) 
PARAMETER (NTB=3, JMX=230, KMX=32, LMX=50) 

DIMENSION XX (JMX, KMX, LMX) , YY (JMX,KMX, IMX) , ZZ (JMX, KMX, LMX) , 

1 R( jinx, kmx, lmx) ,RU(jnx,kinx,lmx) ,RV(jmx,)anx,liBX) , 

2 RW( jmx, kmx, lmx) ,E ( jmx,)anx,lnix) 


WRITE (j to) JMAX, KMAX, LMAX 

WRITE (jto) ( ( ( XX(J,K,L),J« 1,JMAX),K 

* ((( YY (J,K, L) , J= 1,JMAX),K=1 

* ((( ZZ(J,K,L),J= 1,JMAX),K=1 
WRITE (jto) ( ( ( R(J,K,L),J= 1, JMAX) ,K= 

* ( ( (RU(J,K,L) ,J= 1,JMAX) ,K=1, 

* ( ( (RV(J,K,L) , J= 1,JMAX) ,K=1, 

* (((RW(J,K,L) ,J= 1,JMAX) ,K=1, 

* ((( E(J,K,L),J= 1,JMAX),K=1, 


-1, KMAX) ,L*1, LMAX) , 
, -KMAX) , L-’ . LMAX), 

, KMAX) ,L=1, LMAX) 

1, KMAX) ,L=1, LMAX), 
KMAX) ,L-1, LMAX), 
KMAX) , L=1 , LMAX) , 
KMAX) ,L=1, IMAX) , 
KMAX),L=1, LMAX) 
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RETURN 

END 

SUBROUTINE READRES 

PARAMETER fNTB=3 , JMX = 230, KMX— 32, IMX— 50) 

DIMENSION R(JMX,KMX,LMX) ' RU ( 

1 RW(JMX,KMX, LMX) ,E(JMX,KMX,IMX) ,X(JMX,KMX,IMX) , Y ( JKX, KMX, LMX) , 
2 COMMO^DMA/ I NBLK, JMAX(NTB) , KMAX(NTB), LMAX(NTB) 

read restart tape 

NBLK=0 

NTP=8 

50 READ (NTP^ END=99 ) JMAX (NBLK+1) , KMAX (NBLK+1) ,IMAX(NBLK+1) 
NBLK=NBLK+ 1 

WRITE (* ,*) ' READING BLOCK 9 , NBLK 

READ(NTP) ((( X(J,K,L),J= 1, JMAX (NBLK) ) ,K=1,KMAX( NBLK) ) , 

* L 7'( 1 , JMAX (NBLK) ) , K=1 , KMAX (NBLK) ) , 

* **(( 1, JMAX ( NBLK) ) ,K=1, KMAX (NBLK)) , 

* L=1 , LMAX (NBLK) ) 

TRANSLATE AND SCALE ... 

DO 10 J=l, JMAX (NBLK) 

DO 10 K=l, KMAX (NBLK) 

DO 10 L=l, LMAX (NBLK) 

X(J ,K, L) =X( J , K, L) *184 . 

Y(J,K,L)=Y(J,K,L)*184. 

Z ( J,K, L) =Z ( J, K, L) *184 . 

X(J,K,L)=X(J,K,L)+1595.3 
Z (J,K, L) ~Z ( J , K, L) +454.44 
10 CONTINUE 

WRITE (9+NBLK) ((( X(J,K,L) ,J= 1, JMAX (NBLK) ) ,K=1 , KMAX (NBLK) ) , 

* L=l, LMAX (NBLK) ) , 

* ((( Y ( J,K, L) , J= 1, JMAX (NBLK) ) , K=1 , KMAX (NBLK) ) , 

* L=1 , LMAX (NBLK) ) , 

* ((( Z (J,K,L) , J= 1, JMAX(NBLK) ) , K=1 , KMAX (NBLK) ) , 

* L=l, LMAX (NBLK) ) 

READ(NTP)((( R(J,K,L),J= 1, JMAX (NBLK) ) ,K=1, KMAX (NBLK) ) , 

* L=l, LMAX (NBLK) ) , , 

* ( ( (RU(J,K,L) , J= 1, JMAX (NBLK) ) , K=1 , KMAX (NBLK) ) , 

* L=l, LMAX (NBLK) ) , i 

* ( ( (RV(J,K,L) , J= 1, JMAX (NBLK) ) , K=1 , KMAX (NBLK) ) , 

* L=1 , LMAX (NBLK) ) , 

* ( ( (RW(J,K,L) , J= 1, JMAX (NBLK) ) ,K=1, KMAX (NBLK) ) , 

* L=l, LMAX (NBLK) ) , 

* ((( E(J,K,L),J= 1, JMAX (NBLK) ) ,K=1, KMAX (NBLK) ) , 

WRITE (2 9+NBLK) (( ( R(J,K,L) ,J= 1, JMAX (NBLK) ) ,K=1, KMAX (NBLK) ) , 
L=l, LMAX (NBLK) ) , 

( ( (RU(J,K,L) , J= 1, JMAX (NBLK) ) , K=1 , KMAX (NBLK) ) , 

L=1 , LMAX (NBLK) ) , 

( ( (RV(J,K,L) , J= 1, JMAX (NBLK) ) , K=1 , KMAX (NBLK) ) , 

L=1 , LMAX (NBLK) ) , 

( ( (RW (J, K, L) ,J= 1, JMAX (NBLK) ) , K=1 , KMAX (NBLK) ) , 

L=l, LMAX (NBLK) ) , 

( ( ( E(J,K,L),J= 1, JMAX (NBLK) ) ,K=1, KMAX (NBLK) ) , 
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* L=1 , LMAX (NBLK) ) 

GO TO 50 

99 CONTINUE 

DO 5 1=1, NBLK 
REWIND 9+1 
REWIND 29+1 
5 CONTINUE 
RETURN 
END 

SUBROUTINE SCNGEN 

PARAMETER (NTB=3, JMX=230, KMX=32, IMX=50) 
c adjustable : max print, max B.C.'s, max patches per block 

PARAMETER (MPS=10, MBC=25, MP=50) 

c block specific info 
COMMON /CPOIN/ 

* GAMMA , GAMI , GM1R , GGM1 , GSGM , RPR , RPRT , COFMIX , DTBLK , DIS2 , DIS4 , 

* RTPTS, 

* PRESSJ (MBC) , TEMPO (MBC) , 

* PRESSK(MBC) ,TEMPK(MBC) , 

* PRESSL(MBC) ,TEMPL(MBC) , 

* IPWX, IFWY, IPWZ , IPWR, IPWRU, IPWRV, IPWRW, IPWE, IFWQ, 

* IPWS1, IPWVDT, IPWXX, IPWXY , IPWXZ , IPWYX, 

* IPWYY, IPWYZ, IPWZX, IPWZY, IPWZZ, IPWTMU, IPWSPT, IPWCOF, 

* IPW28 , IPW29 , IPW30 , IPWRK1 , IFWRK2 , IPWEPS , IPWAK, 

* NXYZ , IMODE , LTURB , 

* JMAX, KMAX , LMAX, JM, KM, LM, INVISC (3 ) ,LAMIN(3) , 

* NPSEG,JKLPI(3,3,MPS) ,IPORD(2,MPS) , 

* NJSEG , JUNE (MBC) , JKLOW (MBC) , JKHIGH (MBC ) ,JLLOW(MBC) , JLHIGH (MBC) , 

* JTYPE (MBC) , JSIGN (MBC) , INTERJ (MBC) , 

* NKSEG, KLINE (MBC) , KJLOW(MBC) , KJHIGH(MBC) ,KLLOW(MBC) , KLHIGH(MBC) , 

* KTYPE (MBC) , KSIGN (MBC) , INTERK (MBC) , 

* NLSEG , LLINE (MBC) , LJLOW(MBC) , LTHIGH(MBC) , LKLOW(MBC) , LKHIGH (MBC) , 

* LTYPE (MBC) , LSIGN (MBC) , INTERL (MBC) , 

* JEDGE (MBC) , KEDGE(MBC) , LEDGE (MBC) ,JDIR(MBC) ,KDIR(MBC) , LDIR(MBC) , 

* NORMJ (MBC) ,NORMK(MBC) , NORML(MBC) , 

* NJPAT, JPJ2 (MP) , JPJM(MP) ,JPK2(MP) ,JPKM(MP) ,JPL2 (MP) ,JPIM(MP) , 

* NKPAT , KPJ2 (MP ) ,KPJM(MP) ,KPK2(MP) , KPKM(MP) ,KPL2(MP) , KPIM(MP) , 

* NLPAT , LPJ2 (MP) , LPJM (MP) , LPK2 (MP) , LPKM (MP) , LPL2 (MP) , LPLM (MP) 
c 

COMMON /VI S FAC/ ALPHA, BETA,XMACH,C2B, C2BP,VRAT,RE 
NAMELIST/INPUTS/ PREF , TREFR , VRAT , TSUTH , RE , PR, PRT , DDUMP , 

* DTCAF , PCQMAX , SPLEND , SMOO , ST0PL2 , STOPTR , ALPHA , BETA , XMACH , 

* NBLOCK , NMAX , NC , NS PRT , NP , IFXPRT , IFXPLT , L2PL0T , 

* IPLOT , NUMDT , IVARDT , ISOLVE , IRHS , IFILTR, IMUTUR, 

* MBORD , IBORD , NSECTR , LPRES , LREST , GAMMA , DIS2 , DIS4 , 

* NOBORT , IFMAX , IMASS , POsd , TOsd , Lrec4 5 , 

* ORDER, NTURB, IFILTR, NSKIP, STOPTM, CBAR, SSPAN, SREF, 

* REFMX, REFMY, REFMZ 

NAMELIST/BOUNDS/ NJSEG , JUNE , JKLOW , JKHIGH , JLLOW , JLHIGH , 

* NKSEG , KUNE , KJLOW , KJHIGH , KLLOW , KLHIGH , 

* NLSEG , LUNE , ULOW , LJHIGH , LKLOW, LKHIGH , 

i JTYPE, JSIGN, PRESS J,TEMPJ, INTERJ, 

* KTYPE, KSIGN, PRESSK,TEMPK, INTERK, 

* LTYPE, LSIGN, PRESSL,TEMPL, INTERL, 

* JEDGE, KEDGE, LEDGE, JDIR, KDIR, LDIR, 

* NORMJ, NORMK,NORML 
c 

DIMENSION R (JMX, KMX, LMX) ,RU(JMX,KMX,LMX) , RV ( JMX , KMX , LMX ) , 

1 RW (JMX, KMX, LMX) , E ( JMX, KMX, LMX) , X( JMX, KMX, IMX) , Y (JMX, KMX, LMX) , 
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2 Z(JMX,KMX,LMX) 

COMMON/DATA/ NBLK, JMM(NTB) , KMM(NTB) , LMM(NTB) 
LOGICAL FORM 
DATA FORM/. TRUE./ 

DATA MINUS 1/-1/ , IZERO/O/ , IONE/1/ 

BODSYM=0 
INCENT=0 
NOO 
NDUM=0 
NTI = 7 
NTP = 9 
REWIND NTP 
NPNL=0 
ITTR=1 


************************ 


C BASIC DATA 

C** ************** *** *********************** ************* 


READ (NTI, INPUTS) ! Cray version of NAMELIST 

c READ (NTI, NML=INPUTS) ! IRIS version of NAMELIST 

C******** ************************************ *********** 

C PATCH GEOMETRY 

C* ******* **************************************** ******* 


C BODY SURFACE: L=l? J=l, JMX; K-1,KMX 
C SIDE WALL UPPER: K=2; J=l, JMX? L=1,LMX 
r* SIDE WALL LOWER: K=KMX-1; J=l, JMX? L=1,IMX 
END WALL: J=JMX? K=2,KMX-1? L=1,LMX 

CALL BGOM(NP, NFNL, IORDR,Nl, N2 ,N3 , N4 ,N5) 
WHERE: 

NP= PATCH NO. 

NPNL= RUNNING PANEL COUNT 
IORDR= J, K OR L CONSTANT (1,2 OR 3) 
N1 THRU N5= LOOP COUNTERS 
N1,N2 RUN FIRST 
N3,N4 RUN SECOND 
N5 IS CONSTANT COUNTER 


C ASSUMES C-C GRID... 

WRITE (*, *) ' BODY GEOMETRY . . . ' 

NPCH=0 

DO 25 NN=1 , NBLK 

READ (NTI, BOUNDS) ! Cray version of NAMELIST 

c READ (NTI , NML=BOUNDS ) l IRIS version of NAMELIST 

DO 26 J=1,NJSEG 

IF(JTYPE(J) .EQ.7.0R. JTYPE(J) • EQ. 8) THEN 
NPCH=»NPCH+ 1 

IORDR=l - - 

CALL SGOM (NN, NPCH, IORDR, JLLOW ( J) ,JLHIGH(J) ,JKLOW(J) ,JKHIGH(J) , 
1 JLINE(J) ,NPNL) 

ENDIF 

26 CONTINUE 

DO 27 K=1 , NKSEG 

IF (KTYPE (K) . EQ. 7 . OR.KTYPE (K) . EQ . 8 ) THEN 
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NPCH=NPCH+ 1 
I0RDR=2 

CALL SGOM(NN,NPCH, IORDR, KLLOW(K) , KLHXGH(K) ,KJLOW(K) , KJHIGH(K) , 
1 KLINE (K) ,NPNL) 

ENDIF 

27 CONTINUE 

DO 28 L=l, NLSEG 

IF(LTYPE(L) . EQ . 7 . OR . LTYPE ( L) . EQ . 8 ) THEN 
NPCH=NPCH+1 
IORDR=3 

CALL SGOM(NN,NPCH, IORDR, LJLOW(L) , LTHIGH(L) , LKLOW ( L) , LKHIGH ( L) , 
1 LLINE(L) ,NPNL) 

ENDIF 

28 CONTINUE 
25 CONTINUE 

WRITE ( * , * ) ' SCAN DATA INPUT IN FILE SCAN.DAT' 

RETURN 

701 FORMAT (2A4) 

702 FORMAT (515) 

703 FORMAT (6E12. 5) 

704 FORMAT (315/ ( 6E12 . 5) ) 

705 FORMAT ( 6E12 - 5/2E12 .5,15) 

706 FORMAT (6E12. 5) 

707 FORMAT ( 2 0A4) 

END 

SUBROUTINE SGOM (NB , NPTCH , IORDR , N1 , N2 , N3 , N4 , N5 , NPNL) 

PARAMETER (NTB=3 , JMX=230, KMX=32, IMX«50) 

DIMENSION R(JMX f KMX, LMX) , RU ( JMX, KMX, IMX) ,RV(JMX,KMX, HOC) , 

1 RW (JMX, KMX. LMX) ,E(JMX, KMX, HOC) , X(JMX,KMX, HOC) , Y ( JMX, KMX, HOC) , 

2 Z (JMX, KMX. HOC) 

COMMON/DATA/ NBLK, JMAX(NTB) , KMAX(NTB) , LMAX(NTB) 

LOGICAL FORM 
DATA FORM/. TRUE./ 

NTP * 9 
IZERO=0 
ID=2 

NR=N2-N1 

NON4-N3 

NTOT = (NR+1) * (NC+1) 

READ(9+NB) ((( X(J,K r L),J= 1, JMAX(NB) ) , K=1 , KMAX (NB) ) , 

* L=l, LMAX(NB) ) , 

* ((( Y(J,K,L) f J= 1, JMAX(NB) ) , K*1,KMAX(NB) ) , 

* L=1 , LMAX (NB) ) , 

* ((( Z(J,K,L),J= 1, JMAX(NB) ) ,K=1,KMAX(NB) ) , 

* L=l f LMAX(NB) ) 

REWIND 9+NB 

IF (FORM) THEN 

WRITE (NTP f 702) -3, -4,0,0 

WRITE (NTP, 702) NTOT, 1,1,0 

IF ( IORDR . EQ . 1 ) THEN 

WRITE (NTP, 703) ( (X(N5,K,L) ,Y(N5,K,L) ,Z(N5,K,L) ,NB,N5,K,L, 

1 L=N1,N2), K— N3 , N4 ) 

ELSEIF ( IORDR . EQ . 2 ) THEN 

WRITE (NTP, 703) ( (X(J,N5,L) ,Y(J,N5,L) ,Z(J,N5,L) ,NB,J,N5,L, 

1 L=N1,N2), J=N3,N4) 
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ELSEIF ( IORDR . EQ . 3 ) THEN 

WRITE (NTP, 703) ((X(J,K,N5) , Y (J,K,N5) , Z (J,K,N5) ,NB,J,K,N5, 
J=N1,N2), K=N3,N4) 

ENDIF 
ELSE 

WRITE (NTP) -3,-4, 0,0 

WRITE (NTP) NTOT, 1,1,0 

DO 100 I=N3,N4 
IF ( IORDR . EQ . 1 ) THEN 

WRITE (NTP) ( (X(N5,K,L) ,Y(N5,K,L) ,Z(N5,K,L) ,NB,N5,K,L, 

1 L=N1,N2) , K=N3,N4) 

ELSEIF ( IORDR . EQ . 2 ) THEN 

WRITE (NTP) ( (X(J ,N5,L) ,Y(J,N5,L) ,Z(J,N5,L) ,NB,J,N5,L, 

1 L=N1,N2), J=N3,N4) 

ELSEIF ( IORDR . EQ . 3 ) THEN 

WRITE (NTP) ( (X(J,K,N5) ,Y(J ,K,N5) , Z (J,K,N5) ,NB, J,K,N5, 

1 J=N1,N2), K=N3,N4) 

ENDIF 

100 CONTINUE 
ENDIF 

NPNL=NPNL+NR*NC 

701 FORMAT (2A4) 

702 FORMAT (515) 

c 703 FORMAT (6E12. 5) 

703 FORMAT (3F10. 5, 415) 

RETURN 

END 

SUBROUTINE PLTGEN 

PARAMETER (NTB=3 , JMX=2 3 0 , KMX=3 2 , IMX=5 0 ) 
adjustable : max print, max B.C.'s, max patches per block 

PARAMETER (MPS=10, MBC=25, MP=50) 

c block specific info 
COMMON /CPOIN/ 

* GAMMA , GAMI , GM1R , GGM1 , GSGM , RPR , RPRT , COFMIX , DTBLK, DIS2 , DIS4 , 

* RTPTS, 

* PRESSJ(MBC) , TEMPJ (MBC) , 

* PRESSK(MBC) ,TEMPK(MBC) , 

* PRESSL(MBC) ,TEMPL(MBC) , 

* IPWX, IPWY, IPWZ , IPWR, IPWRU, IPWRV , IPWRW , IPWE, IPWQ, 

* IPWS1,IPWVDT,IPWXX,IPWXY,IPWXZ,IPWYX, 

* IPWYY, IPWYZ , IPWZX, IPWZY, IPWZZ , IPWTMU, IPWSPT, IPWCOF, 

* IPW28 , IPW29 , IPW3 0 , IPWRK1 , IPWRK2 , IPWEPS , IPWAK, 

* NXYZ , IMODE , LTURB , 

* JMAX, KMAX, IMAX, JM, KM, LM, INVISC (3 ) ,LAMIN(3) , 

* NPSEG, JKLPI (3,3 ,MPS) , IPORD(2 ,MPS) , 

* NJSEG , JLINE (MBC) ,JKLOW(MBC) , JKHIGH (MBC) ,JLLOW(MBC) ,JLHIGH(MBC) , 

* JTYPE(MBC) , JSIGN (MBC) , INTERJ (MBC) , 

* NKSEG, KLINE (MBC) ,KJLOW(MBC) , KJHIGH(MBC) ,KLLOW(MBC) , KLHIGH(MBC) , 

* KTYPE (MBC) , KSIGN (MBC) , INTERK (MBC) , 

* NLSEG, LLINE (MBC) ,LJLOW(MBC) , LJHIGH (MBC) ,LKLOW(MBC) , LKHIGH(MBC) , 

* LTYPE(MBC) , LSIGN(MBC) , INTERL(MBC) , 

* JEDGE (MBC) , EEDGE (MBC) , LEDGE (MBC) , JDIR(MBC) , KDIR(MBC) , LDIR(MBC) , 

* NORMJ(MBC) ,NORMK(MBC) ,NORML(MBC) , 

* NJPAT , JPJ2 (MP) , JPJM(MP) , JPK2 (MP) , JPKM(MP) , JPL2 (MP) , JPLM(MP) , 

* NKPAT , KPJ2 (MP) ,KPJM(MP) , KPK2 (MP) ,KPKM(MP) ,KPL2(MP) ,KPLM(MP) , 

* NLPAT , LPJ2 (MP) , LPJM (MP) , LPK2 (MP) , LPKM (MP) , LPL2 (MP) , LPLM (MP) 
c 

COMMON /VISFAC/ ALPHA , BETA , XMACH , C2 B , C2 BP , VRAT , RE 
NAMELIST/ INPUTS / PREF , TREFR , VRAT , TSUTH , RE , PR , PRT , DDUMP , 
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* DTCAP, PCQMAX, SPLEND, SMOO , STOPL2 , STOPTR , ALPHA , BETA, XMACH, 

* NBLOCK , NMAX , NC , NSPRT , NP , IFXPRT , IFXPLT , L 2PLOT , 

* I PLOT , NUMDT , IVARDT , ISOLVE , IRHS , IFILTR, IMuTUk, 

* MBORD , IBORD r NSECTR , LPRES , LREST , GAMMA , DIS2 , DIS4 , 

* NOBORT , XFMAX , IMASS , POsd , TOsd , Lrec4 5 , 

* ORDER, NTURB, IFILTR, HSKIP, STOPTM, CBAR, SSPAN, SREF, 

* REFMX, REFMY, REFMZ 

NAMELIST/BOUNDS / NJSEG , JLINE , JKLOW , JKHIGH , JLLOW , JLHIGH , 

* NKSEG , KLINE , KJLOW , KJHIGH , KLLOW , KLHIGH , 

* NLSEG , LLINE , LJLOW , UHIGH , LKLOW, LKHIGH , 

* JTYPE , JSIGN, PRESSJ , TEMPJ, I MTER J , 

* KTYPE , KSIGN , PRESS K , TEMPK, I NTER K , 

* LTYPE , LSIGN , PRESSL, TEMPL, INTERL , 

* JEDGE , KEDGE , LEDGE , JDIR , KDIR , LDIR , 

* NORMJ , NORMK , NORML 

DIMENSION R( JMX, KMX, LMX) ,RU (JMX,KMX,IMX) , RV (JMX, KMX , LMX) , 

1 RW( JMX, KMX, LMX) ,E (JMX, KMX, LMX) , X(JMX,KMX,IMX) ,Y (JMX, KMX, LMX) , 

2 Z( JMX, KMX, LMX) 

COMMON/DATA/ NBLK, JMM(NTB) , KMM(NTB) , IMM(NTB) 

LOGICAL FORM 
DATA FORM/. TRUE./ 

DATA MINUS1/— 1/ , IZERO/O/ , IONE/1/ 

BODSYM=0 
INCENT=0 
NC=0 
NDUM=0 
NTI = 7 
NTP = 9 
REWIND NTP 
NPNL=0 
ITTR=1 


C BASIC DATA 

C******************************************************* 


READ (NTI, INPUTS) i Cray version of NAMELIST 

c READ (NTI, NML= INPUTS) l IRIS version of NAMELIST 
CALL BSIC 


C PATCH GEOMETRY 


BODY SURFACE: L**l; J=1,JMX; K«1,XMX 
SIDE WALL UPPER: K=2 ; J=l, JMX; L=1,IMX 
SIDE WALL LOWER: K-KMX-1; J=l, JMX; L»1,LMX 
END WALL: J=JMX; K=2, KMX-1; L=1,IMX 

CALL BGOM(NP,NPNL,IORDR,Nl,N2,N3,N4,N5) 
WHERE: 

NP= PATCH NO. 

NPNL= RUNNING PANEL COUNT 
IORDR= J,K OR L CONSTANT (1,2 OR 3) 
N1 THRU N5= LOOP COUNTERS 
N1,N2 RUN FIRST 

N3,N4 RUN SECOND 

N5 IS CONSTANT COUNTER 


68 



IF (FORM) THEN 

WRITE (NTP, 701) ' BGIN ' , ' BODG ' 

ELSE 

WRITE (NTP) 'BGIN' , 'BODG' 

ENDIF 

C ASSUMES C-C GRID... 

WRITE ( * , * ) ' BODY GEOMETRY...' 

NPCH=0 

DO 25 NN=1 , NBLK 

READ (NTI, BOUNDS) ! Cray version of NAMELIST 

c READ (NTI f NML= BOUNDS) ! IRIS version of NAMELIST 

DO 26 J=1,NJSEG 

IF (JTYPE ( J) . EQ. 60.0R. JTYPE (J) .EQ.61)THEN 
NPCH=NPCH+ 1 
IORDR=l 

CALL BGOM (NN, NPCH, IORDR, JLLOW ( J) , JLHIGH(J) ,JKLOW(J) , JKHIGH(J) , 
1 JLINE(J) ,NPNL) 

ENDIF 

26 CONTINUE 

DO 27 K=1,NKSEG 

IF(KTYPE(K) . EQ. 60.0R.KTYPE (K) .EQ. 61) THEN 
NPCH-NPCH+1 
IORDR=2 

CALL BGOM (NN, NPCH, IORDR, KLLOW(K) ,KLHIGH(K) ,KTLOW(K) , KJHIGH(K) , 
1 KLINE (K) ,NPNL) 

ENDIF 

27 CONTINUE 

DO 28 L=1,NLSEG 

IF(LTYPE(L) . EQ . 60 . OR. LTYPE (L) .EQ. 61) THEN 
NPCH=NPCH+1 
IORDR=3 

CALL BGOM (NN, NPCH, IORDR, LTLOW (L) , LJHIGH(L) ,LKEOW(L) , LKHIGH(L) , 
1 LLINE(L) ,NPNL) 

ENDIF 

28 CONTINUE 
25 CONTINUE 

IF (FORM) THEN 

WRITE (NTP, 701) 'END ','GRUP' 

ELSE 

WRITE (NTP) 'END ','GRUP' 

ENDIF 

C IF ITTR.GT.O THEN WRITE AERO AWAY TO TAPE 


C AERODYNAMICS 


IF (FORM) THEN 
WRITE (NTP, 701) 
WRITE (NTP, 701) 
WRITE (NTP, 702) 
ELSE 

WRITE (NTP) 
WRITE (NTP) 


'BGIN' , 'AERO' 

'STRT' , 'ITER' 

NPNL , ITTR , I ZERO , IZERO , IZERO 

'BGIN' , 'AERO' 

'STRT' , 'ITER' 
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WRITE (NTP) NPNL, ITTR, I ZERO, I ZERO, I ZERO 

ENDIF 

C CALL AERO ( IORDR, N1,N2,N3,N4,N5) 

p R E • 

C IORDR= J,K OR L CONS TANT (1,2 OR 3) 

C N1 THRU N5= LOOP COUNTERS 

REWIND NTI 

WRITE ( * , * ) ' AERODYNAMICS. . . 9 

DO 35 NN=1 ,NBLK 

READ (NTI, BOUNDS) i Cray version of NAMELIST 

c READ ( NTI , NML= BOUNDS ) ! IRIS version of NAMELIST 

DO 36 J=1,NJSEG 

IF(JTYPE(J) . EQ . 60 . OR. JTYPE (J) .EQ. 61)THEN 
I0RDR=1 

CALL AERO (NN, IORDR, JLLOW (J) ,JLHIGH(J) ,JKU>W(J) ,JKHIGH(J) , 
1 JLINE ( J) ) 

ENDIF 

36 CONTINUE 

DO 37 K=1 , NKSEG 

IF(KTYPE(K) . EQ. 60 • OR. KTYPE (K) .EQ. 61) THEN 
IORDR=2 

CALL AERO (NN, IORDR, KLLOW(K) ,KLHIGH(K) ,KJL0W(K) , KJHIGH(K) , 
1 KLINE (K)) 

ENDIF 

37 CONTINUE 


DO 38 L=1 , NLSEG 

IF(LTYPE(L) .EQ. 60 .OR.LTYPE(L) .EQ. 61)THEN 
IORDR=3 

CALL AERO (NN, IORDR, LJLOW(L) , LJHIGH (L) ,LKLOW(L) , LKHIGH(L) , 
1 LLINE(L) ) 

ENDIF 

38 CONTINUE 
35 CONTINUE 


IF (FORM) THEN 

WRITE (NTP, 701) 'END ','GRUP' 

ELSE 

WRITE (NTP) 'END ','GRUP' 

ENDIF 


C OFF-BODY DATA SCAN 


C C ALL OFFB ( IORDR, N1 ,N2 ,N3 , N4 , N5) 

C w H JFT FT? * 

C IORDR= J,K OR L CONSTANT (1,2 OR 3) 

C N1 THRU N5= LOOP COUNTERS 

C 


IF (FORM) THEN 

WRITE (NTP, 701) ' BGIN • , ' OFFB ' 

ELSE 

WRITE (NTP) 'BGIN' , / OFFB / 

ENDIF 

REWIND NTI 
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WRITE ( * , * ) ' OFF BODY DATA. . . ' 

DO 45 NN=1,NBLK 

READ (NTI, BOUNDS) ! Cray version of NAMELIST 

READ ( NTI , NML= BOUNDS ) ! IRIS version of NAMELIST 

DO 46 J=1,NJSEG 
IF(JTYPE(J) .EQ. 50) THEN 
IORDR=l 

CALL OFFB(NN, IORDR, JLLOW(J) ,JLHIGH(J) ,JKLOW(J) ,JKHIGH(J) , 
1 JLINE(J)) 

ENDIF 

46 CONTINUE 

DO 47 K=1 , NKSEG 
IF(KTYPE(K) . EQ. 50) THEN 
I0RDR=2 

CALL OFFB(NN, IORDR, KLLOW(K) , KLHIGH(K) ,KJLOW(K) , KJHIGH(K) , 
1 KLINE (K)) 

ENDIF 

47 CONTINUE 


DO 48 L=1 , NLSEG 
IF(LTYPE (L) .EQ. 50) THEN 
IORDR=3 

CALL OFFB(NN, IORDR, LJLOW(L) , LJHIGH(L) ,LKLOW(L) ,LKHIGH(L) , 
1 LLINE(L)) 

ENDIF 

48 CONTINUE 
45 CONTINUE 


r 'ID-SPAN PLANE... 

KMID=KMM(l)/2 

IF(2*KMID.NE.KMM(1) )KMID=KMID+1 

CALL OFFB(l,2,2,LMM(l) ,1,JMM(1) ,KMID) 

KMID=KMM(NBLK)/2 

IF(2*KMID.NE. KMM (NBLK) )KMID=KMID+1 

CALL OFFB(NBLK, 2, 1, LMM(NBLK) , 1, JMM(NBLK) ,KMID) 


IF (FORM) THEN 
WRITE (NTP, 701) 
WRITE (NTP, 701) 
ELSE 

WRITE (NTP) 
WRITE (NTP) 
ENDIF 


/ END 9 , / GRUP / 
'END 9 , 'FILE' 

'END ','GRUP' 

'END 9 , 9 FILE' 


c DO 95 J=1 ,NBLK 

C CLOSE 9+J 

C CLOSE 2 9+J 

C 95 CONTINUE 

WRITE (*/*)' PLOT FILE DATA IN FILE PARC-FMT' 
.. 'RETURN 

701 FORMAT (2A4) 

702 FORMAT (515) 

03 FORMAT (6E12. 5) 

/04 FORMAT (315/ ( 6E12 • 5) ) 

705 FORMAT ( 6E12 . 5/2E12 . 5 , 15) 

706 FORMAT (6E12. 5) 

707 FORMAT ( 2 0A4) 
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END 


SUBROUTINE BSIC 

COMMON /VIS FAC/ ALPHA , BETA , XMACH , C2B, C2BP, VRAT , RE 
LOGICAL FORM 
CHARACTER *80 TITL 
DATA FORM/ . TRUE . / 

NTP = 9 
WRITE ( * , * ) 

WRITE (*,*)' PLEASE INPUT TITLE FOR PLOT FILE. ' 

READ ( * , 80) TITL 
RSYM *=0.0 
RGPR =0.0 
ALDEG = ALPHA 
YAWDEG * BETA 
RMACH * 0.0 
RMINF = XMACH 
CBAR *1.0 
SREF *1.0 
SSPAN *1.0 
REFL =1.0 
SCALX * 1.0 
REFMX =0.0 
REFMY = 0.0 
REFMZ =0.0 

IF (FORM) THEN 

WRITE (NTP, 701) 'BGIN' , 'BSIC' 

WRITE (NTP , 707 ) TITL 

WRITE (NTP ,703 ) RSYM , RGPR , ALDEG f YAWDEG , RMACH , RMINF 
WRITE (NTP, 703 ) CBAR* REFL, SREF *REFL* REFL, SSPAN*REFL, REFMX* SCALX , 
1 REFMY *REFL , REFMZ *REFL 

WRITE (NTP, 703) CBAR, CBAR, CBAR, CBAR, CBAR, CBAR 
WRITE (NTP, 701) 'END ' , 'GRUP' 

ELSE 

WRITE (NTP) 'BGIN', 'BSIC' 

WRITE (NTP) TITL 

WRITE (NTP) RSYM, RGPR, ALDEG, YAWDEG, RMACH, RMINF 

WRITE (NTP) CBAR*REFL, SREF *REFL* REFL , SSPAN* REFL , REFMX*SCALX , 

1 REFMY* REFL, REFMZ *REFL 

WRITE (NTP) CBAR, CBAR, CBAR, CBAR, CBAR, CBAR 
WRITE (NTP) 'END ' , 'GRUP' 

ENDIF 

RETURN 

80 FORMAT (A80) 

701 FORMAT ( 2 A4) 

703 FORMAT (6E12. 5) 

707 FORMAT (A80) 

END 

SUBROUTINE BGOM (NB , NPTCH , IORDR , N1 , N2 , N3 , N4 , N5 , NPNL) 

PARAMETER (NTB*3 , JMX=2 3 0 , KMX*3 2 , IMX=50 ) 

DIMENSION R ( JMX, KMX , LMX) , RU (JMX,FMX N< - J/ ngU RgfJM X , K M X^LMX fr 

1 RW (JMX, KMX, LMX) , E( JMX, KMX, LMX) , X (JMX, KMX, IMX) , Y (JMX, KMX, LMX) , 

2 Z (JMX, KMX, IMX) 

COMMON/DATA/ NBLK, JMAX(NTB) , KMAX(NTB), IMAX(NTB) 

LOGICAL FORM 
DATA FORM/. TRUE./ 

NTP = 9 
IZERO=0 
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ID=2 

NR=N2-N1 

NC=N4-N3 

READ (9+NB) ((( X(J,K,L),J= 1, JMAX(NB) ) , K=1,KMAX(NB) ) , 

* L=1 , LMAX (NB) ) , 

* ((( Y (J, K, L) , J= 1, JMAX(NB) ) ,K=1,KMAX(NB) ) , 

* L=1 , LMAX (NB) ) , 

* ((( Z ( J , K, L) ,J= 1, JMAX(NB) ) ,K=1,KMAX(NB) ) , 

* L=1,LMAX(NB) ) 

REWIND 9+NB 


'STRT' , ' PTCH' 

NPTCH , I D , NR , NC , I ZERO 


IF (FORM) THEN 
WRITE (NTP, 701) 

WRITE (NTP, 702) 

IF ( IORDR. EQ. 1) THEN 

WRITE (NTP, 703) ( (X(N5,K,L) ,Y (N5,K,L) , Z(N5,K,L) ,L-N1,N2) 
1 ,K=N3,N4) 

ELSEIF ( IORDR . EQ . 2 ) THEN 

WRITE (NTP, 703) ( (X ( J ,N5, L) , Y ( J,N5 , L) ,Z(J,N5,L) ,L=N1,N2) 
1 ,J=N3,N4) 

ELSEIF ( IORDR . EQ . 3 ) THEN 

WRITE (NTP, 703) ( (X(J,K,N5) ,Y(J,K,N5) ,Z(J ,K,N5) ,J-N1,N2) 
1 ,K=N3,N4) 


ENDIF 

ELSE 

WRITE (NTP) 'STRT' , 'PTCH 7 

WRITE (NTP) NPTCH , ID , NR, NC , I ZERO 

DO 100 I=N3,N4 
IF ( IORDR . EQ . 1 ) THEN 

WRITE (NTP) ( (X(N5,K,L) ,Y(N5,K,L) ,Z(N5,K,L) ,L=N1,N2) , 
1 K=N3,N4) 

ELSEIF ( IORDR . EQ . 2 ) THEN 

WRITE (NTP) ( (X(J,N5,L) ,Y(J,N5,L) ,Z(J ,N5,L) ,Ir=Nl,N2) , 
1 J=N3,N4) 

ELSEIF ( IORDR . EQ . 3 ) THEN 

WRITE (NTP) ( (X(J,K,N5) ,Y(J,K,N5) ,Z(J ,K,N5) ,J=N1,N2) , 
1 K=N3,N4) 


ENDIF 

100 CONTINUE 
ENDIF 

NPNL=NPNL+NR*NC 

701 FORMAT ( 2 A4 ) 

702 FORMAT (515) 

C 703 FORMAT (6E12. 5) 

703 FORMAT (6E2 0.9) 
RETURN 

END 


SUBROUTINE AERO (NB , IORDR , N1 , N2 , N3 , N4 , N5 ) 

PARAMETER (NTB=3 , JMX=2 3 0 , KMX=3 2 , LMX=50 ) 

COMMON/VTSFAC/ ALPHA, BETA, XMACH, C2B,C2BP,VRAT,RE_ _ 

DIMENSION R(oMX,!3!X7l3C??;RW(raX,K!fit,IiMX) ,RV(JMX,KMX,IMX) , 

1 RW(JMX,KMX,LMX) ,E(JMX,KMX,IMX) , X ( JMX, KMX,IMX) , Y (JMX, KMX , IMX) , 

2 Z (JMX, KMX, IMX) 

COMMON/DATA/ NBLK, JMAX(NTB) , KMAX(NTB) , LMAX(NTB) 

LOGICAL FORM 

DATA FORM/. TRUE./, GAMMA/1.4/, GAMI/0.4/ 


NTP = 9 
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GD = GAMMA*GAMI 
FSMACH = XMACH 

CPS - 1./ ( .5*GAMMA*FSMACH**2) 

READ(29+NB) ((( R(J,K,L),J= 1, JMAX(NB) ) , K=1,KMAX(NB) ) , 

* L=1,LMAX(NB) ) , 

* (((RU(J,K,L) ,J= 1 , JMAX (NB) ) , K=1 , KMAX (NB) ) , 

* I>=1, IMAX (NB) ) , 

* ( ( (RV (J, K, L) , J= 1 , JMAX (NB) ) , K=1 , KMAX (NB) ) , 

* L=1,LMAX(NB) ) , 

* (((RW(J,K,L) ,J= 1,JMAX (NB) ) , K=l, KMAX (NB) ) , 

* L=1,LMAX(NB) ) , 

* ((( E(J,K,L),J= 1 , JMAX (NB) ) , K=1 , KMAX (NB) ) , 

* L=1,LMAX(NB) ) 

REWIND 29+NB 


ZER»0. 

I F ( I ORDR . E Q . 1 ) THE N 
DO 5 K=N3,N4-1 
DO 10 L^Nl , N2-1 
RR = 1 ./R (N5 , K, L) 

Q1 = R(N5 ,K, L) 

Q2 - RU(N5,K,L)*RR 

Q3 = RV (N5 , K, L) *RR 

Q4 = RW (N5 ,K, L) *RR 

Q5 = E(N5,K,L) /GAMMA 

PP = GD* (Q5- . 5* (Q2**2+Q3**2+Q4**2) *Q1) 

CPTMP = (PP-1 * ) *CPS 

QQ=SQRT(Q2*Q2+Q3*Q3+Q4*Q4) 

Q2-Q2/FSMACH 
CC Q3-Q3/FSMACH 

CC Q4*Q4/ FSMACH 

IF (FORM) THEN 

WRITE (NTP, 703 )Q2,Q3,Q4,ZER, CPTMP, QQ,R(N5,K,L),Q5 
ELSE 

WRITE (NTP) 02,03,04, ZER, CPTMP, QQ r R(N5 , K, L) , Q5 
ENDIF 

10 CONTINUE 
5 CONTINUE 

ELSEIF ( IORDR. EQ . 2 ) THEN 
DO 15 J=N3,N4-1 
DO 20 L=N1,N2-1 


RR 

= l./R(J/N5,L) 

Q1 

= R(J,N5,L) 

Q 2 

= RU(J,N5,L)*RR 

Q3 

= RV(J,N5,L)*RR 

Q4 

= RW(J,N5,L) *RR 

Q5 

= E(J,N5,L) /GAMMA 

- PP 

= GD* (Q5-. 5* (Q2**2+Q3**2+Q4**2) *Q1) 

CPTMP 

= (PP-1.)*CPS 


QQ-SQRT ( Q2 *Q2+Q3 *Q3 +Q4 *Q4 ) 
Q2~Q 2 /FSMACH 
Cc Q3=Q3/ FSMACH 

CC Q4=Q4/ FSMACH 

IF (FORM) THEN 
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WRITE (NTP, 703) Q2 ,Q3 f Q4 , ZER, CPTMP,QQ,R( J ,N5,L) ,Q5 
ELSE 

WRITE (NTP) Q2,Q3,Q4, ZER, CPTMP, QQ,R(J, N5 , L) ,Q5 
ENDIF 

20 CONTINUE 
15 CONTINUE 

ELSEIF ( IORDR.EQ . 3 ) THEN 
DO 25 K=N3,N4-1 
DO 30 J=N1,N2-1 

= 1 . /R( J , K,N5) 

= R(J,K,N5) 

= RU ( J , K , N5 ) *RR 
= RV(J,K,N5)*RR 
= RW(J,K,N5)*RR 
= E(J,K,N5) /GAMMA 

= GD* (Q5-. 5* (Q2**2+Q3**2+Q4**2) *Q1) 

= (PP-1.)*CPS 
QQ=SQRT(Q2 *Q2+Q3 *Q3+Q4 *Q4 ) 

CC Q2=Q2/FSMACH 

CC Q3=Q3/FSMACH 

CC Q4 =Q4 /FSMACH 

IF (FORM) THEN 

WRITE (NTP, 703 )Q2,Q3,Q4, ZER, CPTMP, QQ,R(J,K,N5) ,Q5 
ELSE 

WRITE (NTP) Q2,Q3,Q4, ZER, CPTMP, QQ,R(J,K,N5) ,Q5 
ENDIF 

30 CONTINUE 
25 CONTINUE 
ENDIF 

703 FORMAT (6E12. 5) 

RETURN 

END 

SUBROUTINE OFFB (NB , IORDR , N1 , N2 , N3 , N4 , N5 ) 

PARAMETER (NTB=3 f JMX=230 ,KMX=32 , LMX=50) 

COMMON/ VI S FAC/ ALPHA , BETA , XMACH ,C2B,C2BP, VRAT , RE 
DIMENSION R(JMX,KMX,LMX) ,RU( JMX, KMX, LMX) ,RV(JMX,KMX,I21X) , 

1 RW(JMX,KMX, LMX) ,E (JMX, KMX, LMX) , X ( JMX,KMX, IMX) ,Y (JMX, KMX, LMX) , 

2 Z (JMX, KMX, IMX) 

COMMON/DATA/ NBLK, JMAX(NTB) , KMAX(NTB) , IMAX(NTB) 

LOGICAL FORM 

DATA FORM/. TRUE./* GAMMA/1.4/, GAMI/0.4/ 

FSMACH = XMACH 
GD=GAMMA*GAMI 

CPS = 1 . / ( 0 . 5*GAMMA*FSMACH*FSMACH) 

NTP = 9 
ZER=0 . 

IONE=l 
IZERO*0 
JN=N2-N1+1 
KN=N4— N3+1 

READ(9+NB) ((( X(J,K,L),J= 1, JMAX(NB) ) ,K=1,KMAX(NB) ) , 

* L=1 , LMAX (NB) ) , 


RR 

Q1 

Q2 

Q3 

Q4 

Q5 

PP 

CPTMP 
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* ((( Y(J,K,L),J= 1 , JMAX (NB) ) , K=l, KMAX(NB) ) , 

* L=1,IMAX(NB) ) , 

* ((( Z(J,K,L),J= 1 , JMAX (NB) ) , K=1 , KMAX (NB) ) / 

* L=1,LMAX(NB) ) 

REWIND 9+NB 


READ(29+NB) ((( R(J,K,L) ,J= 1 , JMAX (NB) ) , K=1 , KMAX (NB) ) , 

* Ir=l,IHAX(NB) ) , 

* ( ( (RU(J,K,L) ,J= 1 , JMAX (NB) ) , K=1 , KMAX (NB) ) , 

* L=1,LMAX(NB) ) , 

* ( ( (RV(J,K,L) ,J= 1, JMAX(NB) ) , K=l, KMAX(NB) ) , 

* ]>1,IMAX(NB) ) , 

* ( ( (RW(J,K,L) ,J= 1, JMAX(NB) ) , K=l, KMAX(NB) ) , 

* L*=1,LMAX(NB) ) , 

* ((( E(J,K,L),J= 1, JMAX(NB) ) , K—l, KMAX(NB) ) , 

* L=1,LMAX(NB) ) 

REWIND 2 9+NB 


IF (FORM) THEN 
WRITE (NTP, 701) 
WRITE (NTP, 702) 
WRITE (NTP, 702) 
ELSE 

WRITE (NTP) 
WRITE (NTP) 
WRITE (NTP) 
ENDIF 


'STRT' , 'BLCK' 

I ONE , IZERO , I ZERO , IZERO , IZERO 
KN , JN , IZERO , IZERO , IZERO 

'STRT' , 'BLCK' 

TONE , IZERO , IZERO , IZERO , IZERO 
KN, JN, IZERO , IZERO , IZERO 


IF (IORDR. EQ.l) THEN 
DO 5 K*N3 , N4 
DO 10 J>N1,N2 

IF(ABS(R(N5,K,L) ) .GT . 1. E-10) THEN 
RR = l./R(N5,K,L) 

ELSE 

RR =1. 

ENDIF 

Q1 = R(N5 f K r L) 

Q2 = RU (N5 , K, L) *RR 

Q3 = RV(N5,K,L) *RR 

Q4 = RW(N5,K,L)*RR 

Q5 = E(N5,K,L) /GAMMA 

PP = GD* (Q5-. 5* (Q2**2+Q3**2+Q4**2) *Q1) 

CPTMP ~ (PP-1.)*CPS 

QQ=SQRT (Q2 *Q2+Q3 *Q3+Q4 *Q4 ) 


IF (FORM) THEN 

WRITE (NTP, 703 )X(N5,K,L) ,Y(N5,K,L) ,Z(N5,K,L) ,Q2,Q3,Q4, 
1 CPTMP, Q5 i 
ELSE 

WRITE (NTP) X(N5,K,L) ,Y(N5,K,L) ,Z(N5,K,L) ,Q2,Q3,Q4, 

1 CPTMP, Q5 
ENDIF 


10 CONTINUE 
5 CONTINUE 


ELSEIF ( IORDR . EQ . 2 ) THEN 
DO 15 J=N3,N4 
DO 20 L»N1,N2 

IF (ABS (R(J,N5,L) ) .GT. 1. E-10) THEN 
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RR 

ELSE 

RR 

ENDIF 

Q1 

Q2 

Q3 

Q4 

Q5 

PP 

CPTMP 


= 1*/R(J,N5 , L) 

= 1 . 

« R( J , N5 , L) 

= RU(J,N5,L)*RR 
= RV(J,N5,L)*RR 
= RW ( J , N5 , L) *RR 
= E(J,N5,L) /GAMMA 

= GD* (Q5- . 5* (Q2**2+Q3**2+Q4**2 ) *Q1) 

= (PP-1.)*CPS 


QQ=SQRT(Q2*Q2+Q3*Q3+Q4*Q4) 

IF (FORM) THEN 

WRITE(NTP,703)X(J,N5,L) ,Y(J,N5,L) ,Z(J,N5,L) ,02,03,04, 
1 CPTMP, Q5 
ELSE 

WRITE (NTP)X(J,N5,L) ,Y(J,N5,L) ,Z(J,N5,L) ,Q2,Q3,Q4, 

1 CPTMP, Q5 
ENDIF 

20 CONTINUE 
15 CONTINUE 

ELSEIF ( IORDR . EQ . 3 ) THEN 
DO 25 K=N3,N4 
DO 30 J=N1,N2 

IF (ABS (R(J, K, N5) ) .GT. l.E-10) THEN 
RR 
ELSE 
RR 


= l.'/R(J,K,N5) 
= 1 . 


ENDIF 

Q1 = R(J ,K,N5) 

Q2 = RU (J, K,N5) *RR 

Q3 = RV(J,K,N5)*RR 

Q4 = RW(J,K,N5) *RR 

Q5 = E(J,K,N5) /GAMMA 

pp = GD* (Q5- . 5* (Q2**2+Q3**2+Q4**2) *Q1) 

CPTMP = (PP-1.)*CPS 

QQ=SQRT ( Q2 *Q2 +Q3 *Q3+Q4 *Q4 ) 

IF (FORM) THEN 

WRITE (NTP, 703) X(J,K,N5) ,Y(J,K,N5) ,Z(J,K,N5) ,Q2,Q3,Q4, 
1 CPTMP, Q5 
ELSE 

WRITE (NTP) X(J,K,N5) ,Y(J,K,N5) ,Z(J,K,N5) ,Q2,Q3,Q4, 

1 CPTMP, Q5 
ENDIF 

30 CONTINUE 
25 CONTINUE 
ENDIF 

701 FORMAT ( 2 A4) 

702 FORMAT (515) 

03 FORMAT (6E 12. 5) 

/03 FORMAT (6E2 0.9) 

704 FORMAT ( 315/ (6E12. 5) ) 

705 FORMAT ( 6E12 . 5/2E12 .5,15) 

706 FORMAT (6E12. 5) 
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707 F ORMA T ( 2 0A4) 
RETURN 
END 



APPENDIX C 


Included here is an example NPARC input file utilized in the zonal aerodynamic analysis of the 
Boeing 727-200 with horizontal and vertical tail leading-edge ice. This example is for the un-iced t-tail 
at zero degrees angle of attack and sideslip, Mach 0.12 and a Reynolds number 1.2 million. Note the use 
of the type 8 boundary condition and the LDIR variable in conjunction with the type 70 boundary 
conditions. The grid structure for this particular problem is summarized below: 


Block 1 

JMAX = 194 

(Freestream direction) 


KMAX = 25 

(Spanwise direction) 


LMAX = 50 

(Normal direction) 

Block 2 

JMAX = 194 

(Freestream direction) 


KMAX = 25 

(Spanwise direction) 


LMAX = 49 

(Normal direction) 

Block 3 

JMAX = 194 

(Freestream direction) 


KMAX = 25 

(Spanwise direction) 


LMAX = 49 

(Normal direction) 


For Block 1, L= 1 coincides with L=2 of Block 3 for so-called wake-cut treatment. For Block 2, K= 1 
and K=KMAX coincides K=KMAX-1 of Block 1 and K= 2 of Block 3 respectively. These overlap 
regions were added automatically with the H-O topology option in PRPAMI. The horizontal and vertical 
tail surfaces are defined by the following indicia! ranges: 


SURFACE 

BLOCK 

J Range 

K Range 

L Range 

Upper vertical 

1 

20 - 170 

1 - 1 

2-35 

Upper horizontal 

1 

64-163 

1-25 

2- 2 

Horizontal tip 

2 

64 - 163 

2-24 

1 - 1 

Lower horizontal 

3 

64-163 

3-25 

1 - 1 

Lower vertical 

3 

20 - 170 

25-25 

1 -49 
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$ INPUTS 

NMAX-5000, NP-3000, 

PREF-15.0, TREFR-500., 

IFXPRT-1, NBLOCK-3, 

DIS2-0.16, DIS4-0.32, IMUTUR = 3, IFTL.TR = 2, 

DTCAP-2.0, PCQMAX-10.0, IVARDT-2, 

NSPRT-50, ST0PL2=l.E-20, ST0PTR-50, STOPTM »' 00002500', 

RE = 1 . 2E+06 , XMACH - 0.12, NOBORT = 1, 

ALPHA - 0.0, BETA = 0.0, 

$END 

SB LOCK :fl 
NPSEG-0 , 

INVISC(1)-1, INVISC(2)=1, INVTSC(3)-1, 

LAMIN(l) -0, LAMIN(2) -0, LAMIN(3) -1, 

SEND 
$ BOUNDS 


NJSEG-2, 


JLINE(1)-1, 

JTTPE(l)-8, 

JSIGN(1)-1, 

JKLOW(l)-2, 

JKHIGH(l)-25, 


JLLOV(l)-l, 

JLHIGH(1)=50, 


PRESSJ(1)-1. 

0, TEMPJ(1)-1.0, 

JTYPE(2)-8, 


JLINE(2)-194, 

JSIGN(2)— 1 

JKLOV(2)-2, 

JKHIGB(2)-25, 


JLLOV(2)-l, 

JLHIGH(2)-50, 


PRESSJ(2)=1. 

0, TEMPJ(2)=1.0, 


NKSEG-5, 

KLINE(1)-1, 

KTYPE(l)-60, 

KSIGN(1)=1, 

KJL0W(l)-20, 

KJHIGH(1)-170, 

KEDGE(l)-25, 

KLLOW(l)-2, 

KLHIGH(l)-35, 

KLINE(2)-25, 

KTTPE(2)=70, 

KSIGN(2)~ 1, 

KJLOW(2)-l, 

KJHIGH(2)-193, 


KLLOW<2)«2, 

ICLHIGH(2)-49, 

INTERK(2)-1, 

KLINE(3)-l t 

KTTPE(3)-50, 

KSIGN(3)-1, 

KJLOW(3)-171,KJHIGH(3)-194, 


KLLOV(3)-l, 

KLHIGH(3)-49, 


KLINB(4)«1, 

KTTPB(4)-50, 

KSIGN(4)-1, 

KJL0V(4)-1, 

KJHIGH(4)=19, 


KLLOV(4)«l, 

KLHIGH(4)-49 t 


KLIH£(5)sl, 

KTTPE(5)=50, 

KSIGN(5)«1, 


KJL0W(5)-20, KJHIGB(5)-170, 

KLL0V(5)-36, KLHIGH(5>-49, 

NLSEG-4, 

LLINB(l)-2, LTTPE(l)-60, LSIGN(1)-1, 

LJL0V(l)-64, UHIGH(1)-163, 

LKL0W(1)-1, LKHIGH(l)-25, LEDGE(l)-45, 

LLINE(2)-1, LTTPE(2)-70, LSIGN(2)-1, 

LJL0V(2)-164, UHIGH(2)«193, 

LKLOW(2)-2, LKHIGH(2)«25, DJTERL(2)-2, LDIR(2)-2, 

LLINE(3)-50, LTTPE(3)«8, LSIGN(3)— 1, 

LJLQW(3)-1, UHIGH(3)=194, 

LKL0W(3)-1, LRHIGH(3)-25, 

PRESSL(3)=1.0, TEMPL(3)=1.0, 

LLINB(4)-1, LTYPE(4)«70, LSIGN(4)-1, 

LJLC7(4)-2, LJHIGH(4)-64, 

LKL0W(4)-2, LKHIGH(4)=25, INTERL(4)-3, LDIR(4)-2, 

$ 

SBLOCK :#2 
NPSEG-O, 

INTISC(1)-1, INVISC(2)-1, INVISC(3)-1, 

LAMIN(l) -0, LAMIN(2) =0, LAMIN(3) =1, 

SEND 

SBOUNDS 
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NJSEG=2, 


JLINE(1)*1, 

JTYPE(l)-8, 

JSIGN(1)=1, 

JKL0W(1)*2, 

JKHIGH(1)»24, 


JLL0V(1)*2, 

JLHIGH(1)=48, 


PRESSJ(1)«1. 

0, TEMPJ(1)-1.0, 


JLINE(2)=194, 

JTYPE(2)-8, 

JSIGN(2)=-1 , 

JKL0W(2)=1 i 

JKHIGH(2)-25, 


JLL0V(2)*1 , 

JLHIGH(2)=49, 


PRESSJ(2)=1 . 

0, TEMPJ(2)=1.0, 


NKSEG-2, 



KLINE(1)=1, 

KTYPE(1)«70, 

KSIGN(1)-1, 

KJL0V(1)~1 , 

KJHIGH(1)«193, 


KLL0V(1)=1, 

KLHIGH(1)=>48, 

INTERK(1)»1, 

KLINE(2)=25, 

RTYPE(2)-70, 

KSIGN(2)— 1, 

KJL0V(2)*1 , 

KJHIGH(2)=193, 


KLL0V(2)-1, 

KLHIGH(2)=48, 

INTERK(2)«4, 

NLSEG-4, 



LLINE(1)*1, 

LTTPE(1)=60, 

LSIGN(1)*1, 

LJL0V(1)=64, 

UHIGH(1)-163, 


uaov(i)«i, 

LKHIGH(1)=25, 

LEDGE(1)=45, 

LLIHE(2)=1, 

LTYPE(2)=82, 

LSIGN(2)-1, 

LJL0W(2)=164 

, LJHIGH<2)-194, 


LKLOV(2)»2 , 

LKHIGH(2)=24, 


LLINE(3)=49, 

LTTPE(3)*8, 

LSIGN(3)— 1, 

LJL0V(3)*1 , 

LJHIGH(3)»194, 


LKL0W(3)=1 , 

LKHIGH(3)-25, 


PRESSL<3)»1. 

0, TEMPL(3)-1.0, 


LLINE(4)*1, 

LTYPE(4)«82, 

LSIGN(4)=1, 

UL0W(4)-1, 

LJHIGH(4)»64, 


LKL0W(4)«2, 

LKHIGH(4)«24, 


$ BLOCK :#3 



NPSEG-O, 



INVISC(1)-1, INVISC(2)«1, INVISC<3)=1, 

LAMIN(l) -0, LAMIN(2) «0, LAMIN(3) -1, 

SEND 



SBOUNDS 



NJSEG.2, 



JLINE(1)»1, 

JTTPE(l)-8, 

JSIGN(1)-1, 

JKLOV(l).l, 

JKHIGH(l)-24, 


JLLOV(l)=l , 

JLHIGH(l)-48, 


PRESSJ(1)*1. 

0, TEMPJ(1)=1.0, 


JLINE(2)*194, 

JTYPE(2)*8, 

JSIGN<2)— 1, 

JKLOW(2)=tl , 

JKHIGH(2)»24, 


JLL0W(2)-1, 

JLHIGH(2>«49, 


PRESSJ<2)»1.0, TEMPJ(2)=1.0, 


NKSEG-4, 



KLINE(1)*25, 

KTYPE(1)=60, 

KSIGN(1)=-1, 

KJLOW(1)*20, 

KJHIGH(1)-170, 


KLLOtf(l)-l, 

KLHIGH(l)-49, 

KEDGE(l)=25 t 

KLINE(2)«1, 

KTYPE(2)«70, 

KSIGN(2)-1, 

KJLOW(2)-l, 

KJHIGH(2)-193, 


KLLOW(2)*l, 

RLHIGH(2)«48, 

INTERK(2)»4, 

KLIKE(3)*25, 

KTYPE(3)-50, 

KSIGN(3)— 1, 

RJLOV(3)«171,KJHIGH(3)*194, 


KLLOW(3)«l, 

KLHIGH(3)*49 f 


KLINE(4)*25, 

KTTPE(4)-50, 

KSIGN(4)— 1, 

RJL0V(4)«1, 

KJHIGH(4)»19, 


KLL0W(4)d t 

KLHIGH(4)«49, 


NLSEG-4 , 



LLINE(1)=1, 

LTYPE(1)=60, 

LSIGN(1)=1, 

LJLOW(l)=64, 

LJHIGH(1)-163, 


LKLOW(l)»l t 

LKHIGH(1)=25, 

LEDGE(1)=45, 
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LLINE(2)«1, LTYPE(2)«70, 

LSIGN(2)-1, 


LJL0V<2)-164, LJHIGB(2)=193 , 

LKL0V(2)*1, LKHIGH<2)-24, 

INTERL(2)»2, 

LDIR(2)*2 , 

LLINE(3)«49, LTYPE<3)«8, 

LSIGN(3)«- 1, 

LJL0V(3)*1, LJHIGH(3)*194, 

LKL0V(3)-1 , LKHIGH(3)»25 , 
PRESSL(3)«1.0, TEMPL(3)d .0, 
LLINE(4)-1, LTTPE(4)-70, 

LSIGN(4)*1, 


ULOV(4)«2, LJHIGH(4)*64, 

LKL0W(4)=1, LKHIGH(4)-24, 

INTERL(4)=3, 

LDIR(4)«2, 


SEND 
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APPENDIX D 


Included here is a discussion of the specific additions made to the NPARC input stream. The user 
should consult the NPARC users manu al 2 as a reference for this discussion. Additional pages to be 
inserted into the VSAERO Revision E.5 User’s Manual are also included. 
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SUMMARY OF NPARC INPUT MODIFICATIONS 


The parameter STOPTM was added to the namelist INPUTS section, (manual, page 69). This 
variable acts much like STOPTR (manual, page 83), but it represents the specific time the job execution 
should terminate normally. The format is in month, day, hour and minute format (MMDDHHMM). This 
option is useful for dedicated time execution on a CRAY computer. 

The boundary condition type 8 was added to foe namelist BOUNDS section which allows the 
specification of external boundary flow-conditions via the data file PARC.SCN (manual, page S3). This 
option acts much die same as the ffeestream type 7 boundary condition except the flow takes on local 
characteristics provided by the PARC.SCN data file. Further information may be found in section 2.3 
of this report. 

The input variables JDIR, KDIR and LDIR were added to die namelist BOUNDS section to 
remove a built-in directionality in the type 70 boundary conditions (manual, page 42). These parameters 
represent arrays s imil ar to the other input within the BOUNDS group (manual, page 57). The user is 
given control of the directionality of a boundary condition region. The values for JDIR, KDIR and LDIR 
are summarized below: 


Segment type 

JDIR 

KDIR 

LDIR 

Value 

Action 

JUNE 

X 



2 

Reverses K index 


X 



3 

Reverses L index 

KUNE 


X 


1 

Reverses J index 



X 


3 

Reverses L index 

LUNE 



X 

1 

Reverses J index 




X 

2 

Reverses K index 


Any other value for these variables will not modify the directionality of the boundary condition region. 
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4.3.5 OFF-BODY VELOCITY SCAN INPUT 


CARD 27: SCAN BOX. 

Cols Variable Value 

1-5 MOLD 

0 

1 

2 

3 

4 

5 

6 
7 
9 


Description 

Shape of velocity scan volume, see 
Sections 3.9.1— 3.9.7 

Stops the scan 

Skewed box option: allows a single 
point, points along a straight line or 
straight lines within a parallelogram 
or within a parallelepiped. 
Requires CARD 28, etc. 

Allows points along radial lines in a 
cylindrical volume. Requires 
CARD 29, etc. 

x,y,z array. Requires CARD 30 

Single panel-to-panel (center-to- 
center) scan. Requires CARD 31 

Multiple panel-to-panel scan. Re- 
quires CARD 32 

Single panel-to-point scan. Re- 
quires CARD 33 

Multiple panel-to-point scan. Re- 
quires CARD 32 and CARD 34 

X,Y,Z array for zonal NS 
applications. Requires CARD 34.1. 


Frmt 

615 
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CARD 27: SCAN BOX. (Cont’d) 


Cols 

Variable 

Value 


Description 


6-10 

MEET 


Controls the intersection line routine (see 
3.9.1) 




Interpolate 

Moved 2 

Flagged 3 

B.L. 4 



0 

Yes 

Yes 

Yes 

No 



1 

No 

No 

No 

No 



2 

Yes 

No 

Yes 

No 



-1 

No 

Yes 

Yes 

No 



-2 

No 

No 

Yes 

No 



-3 

No 

Yes 

Yes 

Yes 

Note: 1 

Interpolation is 

used along the scan line to 

set velocity vectors on points 


inside panel near-field regions using velocities from outside the region 
and surface intersection point velocities. 


Note: 2 Inside points just adjacent to a surface intersection are moved to the in- 
tersection point and take the surface velocity value. 

Note: 3 Points inside bodies are flagged as such. 

Note: 4 Include boundary layer profile if NVPI > 1. 


11-15 


NEAR 


Controls the near-field routine in the ve- 
locity calculation (see 3.9.1) 

0 Active 

1 Active only for surface panels (wake near 

field off) 

-1 Inactive 


Note : Meet of -3 recommended for MOLD 9 option. 
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CARD 34.1: 


GLOBAL COORDINATE ARRAY. (Include if MOLD =9 on 
CARD 27) 


Cols 

Variable 

Value 

Description 

Fnnt 

1-5 

NP1 


Dimensions of array of points entered 

415 

6-10 

NP2 


on CARD 342. If NP1 < 0 array of 


11-15 

NP3 


points read from separate Unit 31 
data file, NAME.SCI. 


16-20 

IPVPRT 

0 

Suppress print of results 




l 

Print results 


75A 

NAME 


Prefix of scan output data file path 
name and optional scan input data 
file path name left adjusted, (e.g., 
NAME.SCO, NAME.SCI). 



CARD 342: 

(Include NP1 • NP2 

• NP3 points) 

• 

Cols 

Variable Value 

Description 

Fnnt 

1-10 

11-20 

21-30 

X(J,KJ.) 

Y(J,KX) 

Z(JJCL) 

Global coordinates of scan point 
Order of input is (((J=1,NP1),K=1, 
NP2),L=1,NP3) 

3F10. 

31-35 

NB 

Block index number transferred to 
NAME.SCO data file. 


36-40 

JP 

Grid point indicial location (J^vL), 
transferred to NAME.SCO data file. 


41-45 

KP 



46-50 

LP 




Note: Format of optional NAME.SCI input data file is a repeat of CARD 34.1 and the 
CARD 342 set This option is useful for very large data input arrays. 
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APPENDIX E 


ENHANCEMENT OF THE DHC-6-200 TWIN OTTER PANEL MODEL 


The existing VSAERO panel model of the DHC-6-200 Twin Otter was refined and enhanced 
under Task 5.3 of the contract. Initially, the existing VSAERO model was reviewed and compared to 
a geometry description of the aircraft supplied by DeHavilland Canada. Several areas of the existing 
panel model were reworked to comply with the available aircraft documentation. These areas included 
die nose of the fuselage, the wing and the vertical and horizontal tail definition. 

The fuselage was first rebuilt, with a continuous maximum width line established through die 
nose. Based upon the DeHavilland documentation, the wing was revised and the intersection with die 
fuselage and the engine nacelle was determined. Similarly, die vertical and horizontal tail was rebuilt 
and a new intersection was defined. The basic refined model is presented in Fig. E.l. 

Further, the wing lift struts, main landing gear strut and wheel and main gear/strut attach-fitting- 
faring were defined. Also, the nose landing gear strut and wheel were added to the model. All struts 
were modeled in their in-flight, extended positions. This enhanced model is shown in Fig. E.2. 

Finally, an actuator disk propeller model was included for simulation of power effects. The wake 
system attached to this disk provides a model of the propeller swirl. As the swirl must have a specified 
rotation, it was decided to represent the propeller on die starboard side of die aircraft. Since a symmetric 
panel model was previously assumed, the port propeller is represented as the mirror image of the 
starboard propeller, which results in a reverse rotation to that on the actual aircraft. This simplification 
may be removed by explicidv defining both port and starboard sides. The resulting model of the DHC-6- 
200 Twin Otter is presented in Fig. E.3. 

At each stage described here, fine-tuning of the panel model was accomplished by examination 
of the VSAERO flow solution and subsequent repair as necessary. 
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Fig. E.2. Front View of 
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Fig. E.3. Rear View of the Enhanced Panel Model of the DHC-6-200 
Twin Otter With Propeller Simulation. 
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